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COMPUTER PROGRAM DESCRIBING TURBINE AERODYNAMIC REQUIREMENTS, 
APPROXIMATE EXTERNAL BLADE GEOMETRIES, AND COOLANT FLOW 
REQUIREMENTS FOR A TWO-STAGE AXIAL-FLOW TURBINE 
by Keith A. Furgalus, David G. Evans, and Michael R. Vanco 

Lewis Research Center 

SUMMARY 

A computer code was written which enables a turbine designer to simultaneously 
determine turbine aerodynamic requirements, approximate external blade geometries, 
and turbine blade and hub and tip wall coolant flow requirements for a two-stage, air- 
cooled, axial-flow turbine . 

Using a free vortex flow model, turbine hub, mean, and tip velocity diagrams are 
generated, along with preliminary inner-stage gas temperatures and pressures. 

Approximate external blade dimensions for each blade row are determined. These 
are based on design data obtained from a survey of existing high-performance aircraft 
gas turbine stator and rotor blades . The resultant designs are able to utilize reasonable 
air cooling schemes, if cooling is required, and at the same time have good aerodynamic 
performance. 

Blade coolant flows are computed by a simplified method based on an assumed 
blade thermodynamic cooling effectiveness. The heat-sink capacity available from the 
cooling air is made to balance the heat sink required to cool each blade row. The blade 
cooling air is adiabatieally mixed with the main gas stream. The innerstage tempera- 
tures and annulus areas are adjusted accordingly. 

Blade hub and tip walls are assumed to be cooled with a tangentially injected air 
film. The inclusion of this wall coolant in the calculations is optional. 


INTRODUCTION 


Several computer codes are available for computing turbine velocity diagrams 
(refs, land 2). Also, several are available which predict turbine blade coolant flow 



requirements and/or metal temperatures. It was felt that there was a need to integrate 
some of the previous efforts in these areas into one package that would enable a designer 
to simultaneously determine turbine aerodynamic requirements, external blade geome- 
tries, and coolant flow requirements. 

A program was written which meets these objectives for a two-stage, air-cooled, 
axial-flow turbine. The program was originally developed to generate designs for a 
turbojet -powered, Mach 3 transport using methane fuel (ref. 3). Based on considera- 
tions of turbine efficiency, size, weight, stress, and cooling requirements, it was felt 
that two-stage turbines were best suited to this application. Hence, the analysis was 
restricted to two-stage designs. The program can, however, be used in any application 
where a two -stage design is desired. 

A complete discussion of the assumptions, equations, computation procedures, and 
computer program is presented along with a complete sample problem. 

A FORTRAN IV program deck is available from the authors. 


ASSUMPTIONS 

The following assumptions were made in the analysis: 

(1) A two -stage turbine is designed. 

(2) A free vortex distribution of the flow from blade hub to tip is assumed in gener- 
ating the velocity diagrams. 

(3) The velocity diagrams are generated assuming a constant turbine mean diameter 
D m , a constant main gas stream specific heat C p , and a constant main gas stream 
specific heat ratio y . The blade coolant airflows are assumed to contribute one-half 
of the equivalent specific work output of the turbine main gas stream. One -half of the 
total blade coolant flow is added at the inlet to the turbine in the calculation of the veloc- 
ity diagrams only. The other half is assumed to bypass the turbine for this computation. 

(4) External blade profiles are determined from design curves obtained from a sur- 
vey of existing stator and rotor blade profiles from high-performance aircraft gas tur- 
bine engines. 

(5) Blade coolant airflows are calculated for each row by the simplified method pre- 
sented in reference 4. 

(6) Blade coolant airflows are mixed with the main gas stream at the exit of each 
blade row. Turbine inner-stage annulus areas and temperatures are adjusted to account 
for the mass flow and energy addition effects of the blade cooling air . 

(7) All flow velocities, flow angles, and blade mean diameters are fixed during the 
calculations of cooling flow requirements to assure no change in turbine work output. 
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(8) Hub and tip wall coolant airflows can be computed if desired. The walls are 
cooled by a tangentially injected air film that is assumed to mix with the main gas 
stream only at the turbine exit. It is also assumed that the wall cooling airflow does 
not contribute to the work output of the turbine. 


PROGRAM DESCRIPTION 

The following description outlines the equations and computation procedures em- 
ployed by the computer program. It is intended to serve as a user’s guide to enable the 
designer to better understand and more effectively utilize the computer program. 

The program is divided into three main sections. The first is the main program 
(TD), which computes the basic turbine velocity diagrams and preliminary inner-stage 
thermodynamic data (i.e. , main gas stream temperatures and pressures) . The second 
is subroutine SURC, which calculates the approximate external blade profiles and sur- 
face areas. The third is subroutine COOL, which computes the required blade coolant 
airflows, mixes these flows with the main gas stream, and iterates until the correct 
combinations of blade cooling flows, inner-stage gas temperatures, and blade heights 
are found. It also computes the required turbine hub and tip wall coolant flows. A de- 
scription of the calculations involved in each of these programs follows. 

Program TD - Velocity Diagram Calculations 

As was previously mentioned, this program calculates the velocity diagrams and 
preliminary inner-stage temperatures and pressures. The velocity diagrams are based 
on a free vortex distribution of the flow from blade hub to tip between each blade row 
using the equations for simplified radial equilibrium: 


V = Constant 

X 

(1) 

r(V u ) = Constant 

(2) 


(Turbine velocity diagram nomenclature is defined in appendix C and fig. 1(a) .) A con- 
stant mean diameter D M , gas specific heat C p g , and gas specific heat ratio y g are 
assumed throughout the turbine. The velocity diagrams are made to satisfy the turbine 
work and hub reaction values specified in the input data. In order to allow for the addi- 
tional annulus area required to pass the blade cooling air, one-half of the estimated 
blade coolant airflow is added to the turbine inlet weight flow to determine the radius at 
which the hub and tip diagrams are determined in the velocity diagram calculations. 
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This cooling air is assumed to contribute one-half of the equivalent specific work out- 
put of the main gas stream. A detailed description of the calculations for each stage 
follows. 

Second stage . - The design begins at the second-stage rotor exit. The exit total 
temperature and pressure are calculated from the given input: 



Ah t, TOT 
C P,g 


( 3 ) 


P h 


2, IV “ P 1,I 


Ah. 


i- 


t, TOT 


c p, g T i } i ?7 TOTy 


(4) 


Note that Ahj. T q T includes the contribution that the blade coolant air makes to the tur- 
bine work output. This work contribution is discussed later, in the description of sub- 
routine COOL. The blade exit tip speed is computed from the input values of exit tip 
radius and turbine rotative speed: 


U T, 2, IV 


n Nr 


T, 2,IV 
30 


Using the specified exit hub-tip radius ratio, the hub radius is computed: 


r H,2,TV 



r T,2,IV 


(5) 


( 6 ) 


The mean radius is computed: 


r M,2,IV 


(r H + r T ) 
_ 


2,IV 


The mean exit whirl is computed: 


V u,M,2,IV 



(7) 


( 8 ) 


4 



where (V u /V cr ) 
computed: 


M, 2,IV 


is specified in the input data. 


The whirl at the hub and tip are 


(V) 


'u,H,2,IV 


m, 2,rv 
r H, 2, IV 


( 9 ) 


V 


(V) 


u,T,2,rv 


M, 2 , IV 
r T,2,IV 


( 10 ) 


The axial absolute velocity V x 2 jy is computed in subroutine VXR: 


'x,2,IV 


w. 


n(r)rdr 
‘'H 


J2,IV 


A check is made to determine whether the computed value of 


(V V cr), 


2, IV 


( 11 ) 


exceeds the 


input limiting value VXEL. If it does, the calculations are terminated and an error mes- 
sage is printed out indicating that the exit tip diameter must be changed to yield a solu- 
tion. If the calculated value is acceptable, the mean and hub blade speeds are computed: 


irNr 


U 


M, 2 , IV 


M, 2, IV 


30 


( 12 ) 


v _* Nr H,2,IV 

Uh ,2,!V ^ 


(13) 


The exit relative velocity at the hub is then computed: 


W H, 2,IV ~ |/ V x,2,IV + ( V u _U ) 


2 

H, 2,IV 


(14) 


The exit relative velocities at the mean and tip radii are computed by substituting the 
mean and tip values for whirl and blade speed into equation (14) . The absolute exit 
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velocities at the hub, mean, and tip are computed: 


V i,2,IV“ V^ V u,i,2,Iv) 2 + ( V x,2,iv) 2 ( 1E 

where i = H, M, T and indicates a substitution of the appropriate values for the hub, 
mean or tip, respectively. 

The calculations for the second-stage rotor inlet are now started. The input value 
of stage reaction is used to compute the relative velocity ratio 


H,IV 


The inlet total temperature and pressure are computed: 


T > _ T , /(2) Ah t,TOT 

A i,iv l 2, rv + \c 

\ p,s 


p i,iv = 


y Ivp’i 


T/^\Ah. 


1 - (2) % TOT 

v c p,g T i,iv% 


where 


T (D Ah t, TOT 


-fe) 


Ti, 


T/^Ahj 
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The stator total-pressure ratio (Pl/Pl) is specified in the input data. The hub-tip 

■ 1 in 

radius ratio at the rotor inlet is not known. However, the required value for W H ^ jy 

is known from equation (16) since W R 2 jy has already been computed. Hence, the 

value of rotor inlet radius ratio which yields the required W g ^ jy must be found. 

Two estimates are made for the inlet radius ratio, denoted in the program as RHRTI1 

and RHRTI2. A check is made to determine whether these radius ratios yield values 

for (V /V 1 which exceed the input limiting value VXIL. The calculation is not 
x cr 1,IV 

terminated if they do. Instead, the smaller of the two estimated radius ratios at the 

inlet is reduced until a solution exists for which (V /V J < VXIL. If no value for 

x cr l,IV 

radius ratio can be found which satisfies this criterion, the input value of VXIL is in- 
ternally increased, and the calculations are restarted at the turbine exit. Having two 

suitable inlet radius ratio estimates, a curve of Wt, •, n r against (r„/ r™) is gen- 

n, i, iv xi A iiv 

erated by subroutine MAXIM. The curve is generated by computing Wjj ^ jy from 


W- 


H, 1,IV 




vy 


H, 1,IV 


+ v: 


x , 1,IV 


( 21 ) 


where using values for inlet hub-tip radius ratio between RHRTI1 and RHRTI2 and letting 
r M l iv = r M 2 IV’ anc * blade speed are computed: 


V u,H,l,IV 


Ah. 


t, TOT 


+ — (V r) 
30 U 


M, 2 , IV 


^■ r H,l,IV 


( 22 ) 


U H,1,IV 


_ r 

- r H,l,IV^ 


(23) 


The inlet axial absolute velocity is computed in subroutine VXR as before: 


x, 1,IV 


w. 


n 

it I 

J H 


p(r)rdr 


i,rv 


(24) 


Subroutine MAXIM then proceeds to find the minimum value of W H ^ jy and its cor- 
responding radius ratio. This newly computed radius ratio is used as the new minimum 
radius ratio, replacing the previously estimated minimum. If these limiting inlet radius 
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ratios, RHRTI1 and RHRTI2, yield values for W H ^ jy which bracket the desired value 
computed from equation (16), subroutine ROOT proceeds to find the value of radius ratio 
between RHRTI1 and RHRTI2 which yields the desired value for Wjj j jy. If a solution 
is not possible with the previously computed radius ratio limits, the exit tip radius is 
changed and the calculations restarted. Having found the correct inlet radius ratio, the 
following velocities are computed at the rotor inlet: 



1,IV 


= — (r.) 

30 1 1,IV 


(25) 


jtN 


7 (|) Ah t, TOT + ~ <V> 


'u,i,l,IV 


M, 2,1V 


jtN 

30 


(26) 


r. 


i,l,IV 


W. 


i,l,IV 




IV*(V U -C) 




(27) 


V i, 1,IV ~ y< V u,i, 1,IV )2 + (V x, 1,IV^ 


(28) 


First stage . - The calculations are now continued for the first- stage rotor. The 
rotor inlet and exit total pressures and temperatures are computed: 


Po 


2,11 


_ P 1,IV 


(29) 


HP * „ rri t 

A 2,n ~ i i,iv 


T « _ T « . r 0D Ah t,TOT 

L i,n ~ l 2,u 

p,g 


(30) 

( 31 ) 
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(32) 


P 


i,n 


p 2 ,n M\ 

y \Ply/ 

T (D Ah t,TQTY 1 

c p,g T i,n 1? (D/ 


The first-stage velocity calculations are similar to those for the second stage. The 
major difference in the calculations for the first stage is that the mean exit whirl has not 
been specified as it was for the second stage. Two practical design limits between 
which a solution should exist are chosen: 


and 


V u,M,2,n = (°- 1V cr> 


M,2,n 


(33) 


V u,M,2,n = (-°- 4V cr) 


M,2,n 


(34) 


The annulus area at the first-stage exit is initially assumed equal to the area at the 
inlet to the second-stage rotor times the area ratio across the second-stage rotor. The 
first-stage rotor exit radii are then computed: 


r T,2,II ~ r T, l,IV 


y AR rv 


(35) 


r M,2,IX ~ r M,l,IV 


(36) 


r H,2,n ~ 2r M,2,n ' r T,2,n ^ 

Using the first assumed value for mean exit whirl (i . e . , 0.1 V cr M j j), exit blade 
speeds and whirls are computed at the hub, mean, and tip by the same procedure em- 
ployed for the second stage. The axial absolute velocity ratio (V /V" ) is computed 

x cr 2 n 

as before in subroutine VXR. If the computed value exceeds the input limiting value 
VXEI, the area ratio is increased until an acceptable value for the axial velocity ratio is 
found. The exit hub relative velocity is then calculated from 


W 


H,2,n 




(V x 2 n )2 + (V n - U > 2 

u H,2,n 


(38) 
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The calculation proceeds to the inlet where the hub, mean, and tip radii are computed: 


r T,i,n ~ r T,2,n V AR n ^ 


r M,i,n ~ r M,2,n ( 4 °) 


r H , l , n - 2r M, 1 , n " r T , 1 , n (41) 


The blade hub, mean, and tip inlet blade speeds and whirls are computed exactly as 

before. The axial velocity ratio (V„/V„ J is computed in subroutine VXR and 

x cr 1,H 

checked against its input limit VXII. Again, the area ratio is increased, if necessary, 
to yield an acceptable axial velocity ratio . The inlet hub relative velocity Wjj ^ jj is 
then computed as before. With known values for ^ jj and W H ^ jp ^e stage hub 
reaction can be computed: 



(42) 


The calculations are repeated for the second assumed mean exit whirl 
(i.e. , -0.4 V cr ^ 2 jj)' These calculations yield a second value for the stage hub 
reaction. If these two calculated values for hub reaction bracket the desired value sup- 
plied in the input data, a series of calculations similar to those in subroutine ROOT finds 
the value of mean exit whirl between the two assumed design limits which yields the de- 
sired stage hub reaction. If the two initially assumed values for mean exit whirl do not 
yield a solution, the calculations are terminated with an error message indicating that 
in the mean exit whirl iteration there is no solution and that the hub-tip radius ratio or 
whirl limits must be changed. 


Subprogram SURC - External Blade Geometry Calculations 

A method was developed to determine the approximate external blade configurations 
required for each blade row in the turbine. The method consisted of surveying existing 
government and industrial high-performance aircraft gas turbine stator and rotor blade 
profiles. From the survey, a series of working curves were made which relate the re- 
quired design parameters. These curves were then used to obtain the dimensional char- 
acteristics of blading required by interpreting the curves at the mean radius values of 
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blade gas entering and leaving angles 0^ and 0 2 . (Nomenclature is defined in appendix C 
and fig. 1(b) .) These design curves are supplied as an integral part of the program. 

As noted, the starting point in the calculation is that of determining the blade enter- 
ing and leaving angles 0^ and 0 2 at the mean radius from the data generated in pro- 
gram TD. With these angles, an initial blade camber angle <p is found by interpolation 

z 

from figure 2, which is based on data from references 3 and 5. This interpolation is 
performed using a four -point bivariate technique (subroutine BVTL). 

The value for the stator aerodynamic loading coefficient x(/ is interpolated from 
figure 3. The loading coefficient for the rotors is set equal to 1. From the survey, it 
was found that the blade camber angle was influenced by the aerodynamic loading coeffi- 
cient and by the particular row for which the blading was designed. The effect is shown 
in figure 4, which relates i[/ to a correction A (p for each of the four blade rows . The 
final value for the blade camber angle for each blade row is then computed from 


<P = (p z + A(p 


(43) 


Using these values for camber angle along with a A/3 defined as 


A/3 =j3j_ -0 2 


(44) 


values for blade camber length to axial chord length ratio (C T / C ) are found for each 

J-, x z 

row from figure 5; the data being taken from references 3 and 5. Blade chord lengths 
are determined by dividing the blade heights found in program TD by their aspect ratios, 
which were found, from the survey, to equal 2, 3, 3, and 5, respectively, for the four 
blade rows. In cases where these values of aspect ratio result in an axial chord length 
of less than 1 inch (2.54 cm) for any of the blade rows, the aspect ratios for all rows are 
reduced and the chord lengths recomputed for all the rows until the minimum value of 
axial chord length for any row is greater than 1 inch (2 . 54 cm) . This is done to assure 
that the cross-sectional area of the various blade profiles will be large enough to accom- 
modate internal cooling air passages. 

Note that (C T /C ) for the stators is corrected as a function of blade thickness 
L x z 

ratio t^„ /C T , as indicated in figure 6. A value for blade thickness ratio of 0. 16 was 
assumed for blading developed in this program. No corrections are applied to the rotor 
blades . 

It was also noted that the surveyed blading had airfoil perimeter lengths approxi- 
mately 2.03 times their mean camber lengths at blade thickness ratios of 0. 16. Thus, 
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a single blade external surface area Ag b is computed for each blade in a row: 


( 45 ) 


A S b = 2.03 (C L ) (BH) 

The number of blades in each row is found by dividing the mean radial circumference (a 
constant throughout the turbine) by the blade pitch for that row, where the pitch is com- 
puted from 


P = 


*< C x> 

2 sin^ ^ (cot 0g - cot 0j) 


(46) 


With this information and the velocity diagram angles and blade heights generated 
in TD, approximate external blade profiles are described. The resultant external blade 
surface areas are used as the starting point in the iteration to determine blade coolant 
airflow requirements (if cooling is required), as described in the next section. 


Subprogram COOL - Turbine Blade and Wall Coolant Flow Calculations 

The coolant airflow required to cool both the turbine blades and the hub and tip walls 
is computed next . 

Rather than becoming concerned with the complexities associated with individual 
blade cooling configurations and determinations of cooling airflow requirements based 
on local metal temperature levels and heat-transfer coefficients, the simplified method 
described in reference 4 is employed. The method consists of making a heat balance 
between the heat sink required to cool each blade row to the heat sink available in the 
cooling air . The method requires knowledge of the average main gas stream conditions 
in each row, the average blade midspan metal temperatures, the cooling air supply tem- 
perature, and the thermal cooling effectiveness of the blade cooling design. 

To start the coolant calculations, a base design is needed. This design is obtained 
from the previously computed velocity diagrams and blade external surface areas. All 
flow velocities and angles, and the turbine mean diameter computed in the base case are 
fixed during the coolant flow calculations to assure no change in turbine work output. 

The calculations for the base turbine assume that one-half of the initially estimated 
blade coolant flow is added to the burner discharge flow and that this flow remains con- 
stant throughout the turbine and does not affect the turbine inlet temperature . The other 
half of the coolant flow bypasses the turbine. In subsequent cooling iterations, the total 
coolant flow from each blade row is mixed with the main gas stream at the exit of that 
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blade row, and the turbine inner-stage weight flows and temperatures are adjusted 
accordingly . 

Using the results of the base case turbine velocity diagram and blade surface area 
computations, a calculation is performed to determine individual blade coolant require- 
ments. The blade gas-side heat transfer coefficient is computed: 


h = 0 , 037 

O 




(47) 




g 


"L,M 


where the velocity V is the velocity relative to the row. The value p V is the aver- 

& . o o 

age between the row inlet and exit. The transport properties are "based on the average 
total relative temperature of the gas entering the blade row. The required heat sink is 
computed for each blade row: 


Q R = h g A S,b,TOT ( T g 


BM> 


(48) 


Note that the average midspan blade metal temperature Tgj^ has been specified in the 
input data . 

To obtain the heat sink available from the coolant Q^, we use 



(49) 


Modifying this to include the thermal cooling effectiveness of the blade cooling configur- 
ation, 


(Tg - T j) 




y,b 


rn rn TT 

i BM " -4,y,b 


(50) 


we get 


Qa = % < wC p) h ^ T BM 

y ? ° 


- T l y.b> 


(51) 
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Since for each blade row, 

W y,b = \^,b,TOT^ "^bm* (52) 

W E W E (C p* h ( T BM - T i,y,b) 

If any blade row in the turbine does not require cooling, the program will calculate zero 
cooling for that row. 

This procedure assumes no circumferential gas temperature differences. In prac- 
tice, the first- and second-stage stators are exposed to these gradients. This tempera- 
ture variation is accounted for by multiplying the required coolant flow determined in 
equation (52) by the ratio 


rp TT 

g, max 


" T BM 


T - T 

l g 1 BM 


(53) 


where T” max is the maximum circumferential hot-spot gas temperature. For the 
first-stage stators 


T" 

g, max, I 


T" + PF[?'-(T y;1+ AT y § 


(54) 


where T AT , and PF are specified input data, and are respectively the tempera- 
y j y 

ture of the coolant air entering the turbine, the difference between this temperature and 
the compressor discharge temperature, and the combustor pattern factor. For the 
second-stage stators, the maximum hot-spot gas temperature is determined by assuming 
that one-half of the increase in gas temperature in the first-stage stator due to the effect 
of pattern factor is attenuated across the first stage before entering the second stage so 
that 


(T" - T") 

- v g, max g'j j 

m M _ »T» rf . x ? 1 

A g, max, in g,HI + «, 


(55) 


The coolant flows computed by this method provide for cooling of the stator location ex- 
posed to the circumferential hot-spot gas temperature. The remaining stators are 
overcooled. 

After having computed the individual blade row coolant flows, the turbine flow is 
recomputed at each location in the turbine taking into account the addition at each loca- 
tion of the blade coolant flow of each blade row upstream. The coolant flows are 
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adiabatically mixed with the main gas stream at each blade row exit: 


T 


g,2 


(w t T'p i + (wT") yb 


w. 


t,2 


( 56 ) 


These newly computed temperatures are then used to compute new values for p V_ 
across each row. The revised blade surface areas are then computed: 


A S,b,TOT,n+l 


' A S,b,TOT,n 


W, 

(w t , n ) <pV) n+1 


(57) 


(For the rotors, V is replaced by W.) These steps are repeated until successive iter- 
ations yield repeating combinations of blade cooling airflow ratios, inner- stage gas tem- 
peratures, and required external blade surface areas for each row. 

The wall hub-tip cooling airflow ratio w y w /w E is computed next. This calcula- 
tion is optional (see input instructions) . The cooling method employed is that of film 
cooling air injected tangentially to the hub and tip walls upstream of each blade row. 

The following correlation is used to predict wall coolant flow ratios: 


w. 


L914 C^g^O.2 (p g V g C L ) 0 - 8 D m 


y,w 


w 


E 


w 


E 


p,y/ 

1.9 Pr 


0.67 


(58) 


where 


1600 - (T" - 150) 
Ty - (T’g - 150) 


(59) 


The development of this equation is presented in reference 3. It is assumed that the 
wall .cooling air remains as a cool boundary layer along the walls. As a result, the wall 
cooling air is assumed not to reduce the inner-stage main gas stream temperatures or 
to produce any work in subsequent turbine stages. The wall coolant is assumed to mix 
with the main gas stream only at the turbine exit. A final adjustment of blade height in 
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each row is required to accommodate the hub and tip wall coolant airflows. Assuming 
that the injection velocity of the wall coolant is one-half of the absolute axial velocity of 
the main gas stream at the exit to each row and that the coolant static pressure is equal 
to the main gas stream static pressure at the exit of each row, the blade height adjust- 
ment due to wall coolant flow is computed: 


(ABH) 


/ w 

\P% 


y,w 


y,w 


(BH) g+y,b A a,(g+y,b) 


(60) 


This additional length of blading is assumed to be cooled by the wall cooling air. There- 
fore, an increase in blade cooling airflow is not needed. 

The preceding calculations yield a complete design based on the initially assumed 
values of coolant flow ratios for the blades and walls. If the final calculated values for 
these coolant flow ratios do not equal the initially assumed values which are provided 
as input data, the program must be recycled. Since the turbine specific enthalpy change 
is given 



(i . w y,TQT \ 
\ W E / 


(1 + f/a) 



r 

w y,b 

Li 

W E 


2 



\ W E ) 


- 



(61) 


a change in coolant flow ratio affects the turbine work output (assuming that the compres- 
sor enthalpy change or turbine shaftwork stays constant). Since the velocity diagrams are 
computed using the input assumed coolant airflow ratios, the velocities used in the pro- 
gram are not correct if the assumed and calculated coolant flow ratios do not agree. 

Thus the program must be recycled using the newly calculated values as input until the 
calculations predict the same coolant flow ratios as were assumed in the input. A cor- 
rect overall solution is then assured. All the iterations leading to a correct solution are 
printed out by the program. 
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Utility Subprograms 

Several utility subroutines are used by TD, SURC, and COOL. These are briefly 
described below. 

Subroutine ROOT . - Working between two end points a and b, this subroutine finds 
the value of X which yields the desired value of Y for an arbitrary function Y = f (X) 
where a ^ X < b. It is used in the program to find the radius ratio which yields the 
desired hub inlet relative velocity for the rotors. 

Subroutine RHRT2 . - This subroutine calculates the axial absolute velocities at the 
mean radius and compares them to the limiting values prescribed in the input data. It 
reduces the hub-tip radius ratios in an effort to reduce the axial velocities if the calcu- 
lated values exceed the specified limits. In the case of the second-stage rotor inlet, the 
subroutine will increase the input limit VXIL if necessary, and return to the beginning 
of the calculations . 

Subroutine DENR. - The axial velocities V v at the mean rauius are computed based 
on compressible flow. 

Subroutine VXR . - This subroutine calculates densities at the hub, mean, and tip 
and then integrates to find the axial velocities, 

Subroutine MAXIM . - This subroutine is used to find the minimum value on a curve 
of relative hub velocity against radius ratio. 

Subroutine BVTL . - This is an interpolation subroutine used to interpolate values 
when the desired value actually lies on a surface (i.e. , a family of curves are presented 
from which the desired value must be chosen as a function of two other coordinates) . 

The procedure is described in reference 6 under the subject heading "Four Point Bivar- 
iate Interpolation.” An error message - "Bivariate Interp Error VI = XXX. X, 

V2 = XXX. X” indicates that the coordinates presented to the subroutine are out of range 
of the data provided for the particular set of curves. The subroutine does not extrapo- 
late data. 

Subroutine LINT . - This is a standard linear interpolation routine used whenever 
data are required from a curve which is described by y = f(x) . 


DESCRIPTION OF REQUIRED INPUT DATA 

The following instructions provide a detailed description of the required program 
input data. The "Line" designation refers to one data card. The "Location" specifies 
the card columns in which the input variable must be entered. The "Type of Number" 
refers to either integer or fixed point numbers. An integer may consist of more than 
one digit but must be oriented to the right in its field. A decimal point must not appear 
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in the field. A fixed point number may appear anywhere in its field. A decimal point 
should appear in the field if the number is not right adjusted. All the data cards must be 
in the deck. If an uncooled turbine is being considered, the data cards must still contain 
cooling input data; the program will predict zero coolant flow based on the input data. 

An approximate initial input value for DHT can be computed by equation (61) using 
the estimated cooling flow ratios and a known compressor enthalpy change or turbine 
shaft work requirement. 

If wall cooling is to be ignored (i . e. , NWC = 1), the values entered for Y and YB 
should be equal. 


Line Location Type of FORTRAN Description 




number 

symbol 


1 

1-3 

I 

NCASE 

Identification number for the case 

2 

1-12 

FP 

TCI 

Temperature of the coolant air as it enters the 
blades and hub and tip walls, °R 

2 

13-18 

FP 

CORRF 

A value of 1 . 0 must be entered 

2 

19-24 

FP 

CORRFW 

A value of 1 . 0 must be entered 

3 

1-12 

FP 

FAR 

Primary burner fuel-air ratio 

3 

13-24 

FP 

Y 

Initially estimated value of the total (blades 
plus wall) coolant airflow ratio Wy tqt^E 

3 

25-36 

FP 

YB 

Initially estimated value for the blade coolant 
flow ratio w TI K /wp 

4 

1-6 

FP 

TTO 

Turbine total inlet temperature, R 

4 

7-12 

FP 

PTO 

Turbine total inlet pressure, lbf/ft^ abs 

4 

13-18 

FP 

RPM 

Turbine rotational speed, rpm 

4 

19-24 

FP 

DRTE 

Exit tip radius increment, in. (A value of 0. 1 
is usually sufficient.) 

4 

25-30 

FP 

WO 

Turbine inlet gas weight flow, lbm/sec 

4 

31-36 

FP 

DHT 

Turbine enthalpy drop, including the contribu- 
tion from the blade cooling air (see eq. (61)), 
Btu/lbm 

4 

37-42 

FP 

RTE 

Assumed tip radius at the second-stage exit, 
in. 

Assumed hub-tip radius ratio at exit of last 
stage 

4 

43-48 

FP 

RHRTE 

4 

49-50 

I 

NWC 

NWC = 1, if only blade coolant flows are to be 
calculated; NWC = 2, if both blade and wall 
coolant flows are to be calculated 
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Line 

Location 

Type of 
number 

FORTRAN 

symbol 

Description 


5 

1-6 

FP 

WS(1) 

Work fraction for first stage 


5 

7-12 

FP 

WS(2) 

Work fraction for second stage 


6 

1-6 

FP 

GAMMA 

Main gas stream specific heat ratio 

^0* 

6 

7-12 

FP 

R 

Main gas stream gas constant, ft-lbf/(lbm)( 0 R) 

6 

13-18 

FP 

PC 

Stator total -pressure ratio (p^/pp 


7 

1-6 

FP 

ETAT 

Overall turbine total efficiency 


7 

7-12 

FP 

ETA(l) 

Total efficiency of first stage 


8 

1-6 

FP 

R(l) 

Hub reaction for first- stage rotor 




8 7-12 FP R(2) Hub reaction for second-stage rotor 



9 

1-6 

FP 

VUMCR 

Mean exit whirl of last-stage rotor 

CVu/Vcr) „ 
u cr M,2,IV 

10 

1-6 

FP 

VXEL 

Limiting value on axial velocity ratio V x /V cr 
at exit of second- stage rotor 

10 

7-12 

FP 

VXEI 

Limiting value on axial velocity ratio V /V 

x cr 

at exit of first-stage rotor 

10 

13-18 

FP 

vxn 

Limiting value on axial velocity ratio V x /V cr 
at inlet to first- stage rotor 

10 

19-24 

FP 

VXIL 

Limiting value on axial velocity ratio V /V 

A vl 

at inlet to second- stage rotor 

11 

1-12 

FP 

FARC 

Value for FARC must be set equal to value for 
FAR specified in line 3 

11 

13-24 

FP 

PF 

Combustor pattern factor 

11 

25-36 

FP 

DTY 

Decrease (+DTY value) or increase (-DTY 
value) in temperature of cooling air between 
compressor bleed and inlet of blade or wall 
coolant passage 

12 

1-3 

I 

NPTD 

Number of temperatures at which gas stream 


fluid properties are entered 
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Line 

Location 

Type of 
number 

FORTRAN 

symbol 

Description 

13 

1-12 

FP 

TPR 

First temperature at which fluid properties 
will be entered, °R 

13 

13-24 

FP 

PRN 

Value of Prandtl number at first temperature 
entered 

13 

25-36 

FP 

VIS 

Value of viscosity at first temperature entered, 
lbm/(sec)(ft) 

13 

37-48 

FP 

CON 

Value of thermal conductivity at first tempera- 
ture entered, Btu/(sec)(ft)(°R) 


The remaining values of temperature, Prandtl number, viscosity, and thermal conduc- 
tivity are entered on (NPTD-1) cards according to the format just described. 


13+NPTD 

1-12 

FP 

TMS 

Average stator midspan metal temperature, °R 

14+NPTD 

1-12 

FP 

TMR 

Average rotor midspan metal temperature, °R 

15+NPTD 

1-12 

FP 

THERME 

Turbine blade thermal cooling effectiveness 
(This value is used for all four blade rows.) 

16+NPTD 

1-12 

FP 

CPC 

Average specific heat for cooling air, 
Btu/(lbm)(°R) 

17+NPTD 

1-12 

FP 

TWALLA 

Average blade hub and tip wall metal tempera- 
ture, °R 


DESCRIPTION OF COMPUTER OUTPUT DATA 

Several sets of computer output are presented in appendix A . Each set contains 
three sections - ’’velocity diagrams, ’’ "initial surface area calculations, " and "coolant 
calculations. " Each set is the result of one iteration to converge to a correctly matched 
solution, as described previously. This allows the designer to examine intermediate re- 
sults to observe trends in the iteration. The format in all the sets is identical. Hence, 
only one set will be described. 

In the velocity diagrams and coolant calculations sections, a subsection titled 
"input data" is printed out. This is a recapitulation of data read and/or used by routines 
TD and COOL. The FORTRAN symbols and units are the same as those used in the input 
instructions. The units used in the printout are defined in table I. The velocity dia- 
grams section of the printout is the result of the computations in program TD. The inlet 
and exit designations refer to the rotor inlet and exit. Since total and total relative data 
are presented, the stator data is also available. The actual velocity diagrams are 
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described under the headings RADE, BAE, VAE, etc. The first line of data under these 
headings is data at the turbine tip radius. The second and third lines are for the mean 
and hub data, respectively. (The heading nomenclature is related to the velocity dia- 
grams in fig. 7.) The numbers printed directly under the velocity ratios are the actual 
velocities at that location (i.e. , the numerator in the velocity ratio) . Note that when the 
individual blade row cooling air is mixed with the main gas stream later in subroutine 
COOL, different inner-stage temperatures and hence, different critical velocities re- 
sult. Thus, the velocity ratios printed out here will change. The actual velocities, 
however, do not change. The critical velocities computed after the addition of the blade 
cooling air are printed out in the coolant calculations section. Thus, using the actual 
velocities and the revised critical velocities, corrected velocity ratios can be computed. 

The data presented in the initial surface area calculations section are the results 
from subroutine SURC. Blade geometry symbols are defined in figure 1(b). The values 
for blade surface areas are those for the base case which are later updated. 

The data presented in the coolant calculations section are the results from subrou- 
tine COOL. The word "initial” in the descriptions refers to the base case in which one- 
half of the assumed blade coolant flow is added at the turbine inlet and the other half 
bypasses the turbine (i.e. , the assumption used to generate the velocity diagrams). 

The word "revised" refers to values computed after the blade cooling air from each row 
has been mixed with the main gas stream. Note that for total relative temperatures at 
the rotor inlet or exit, the values presented are at the mean radius. All coolant flow 
ratios are obtained by ratioing the cooling air weight flow to the engine inlet weight flow. 
If wall cooling is included (i.e. , NWC = 2), the blade heights and turbine exit tempera- 
ture printed out have been corrected to include the effects of wall cooling airflows . 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 3, 1970, 

720-03. 
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Sample Problem Output 


CASE-144 

VELOCITY DIAGRAMS ********** 
INPUT DATA** 


TC I =1204 . 900 
Y= 0.076' 
RPM-6325. :;00 
R TE = 22.000 

N VUM= 1 
GAMMA= 1.531 
ETAT = 0.903 

VUMVCR= -0.039 
VX I L= 0.550 


CGRRF = I.CoO 

YL- C.C430 
DRTE = 0* ICO 

RHRTE- 0,641 
NRTE = I 
Ct>= 0.300 
ETAI 11- 0.890 

V XE L- 0.8 0 v 


CGRRFw= .1* COO 

TTO=2660.0. 0 
W0= 316. 7c: 

N= 2 

KSCIi* 0.360 
R= 53.960 
RE 111= 0.180 

VXF I = 0.6GC 


FAR= C.023 
■PTO-18734.- 30 
OH 7 = 141.900 
NR— 2 

W5t?)= 0.44i. 

PC= 0.970 
RE { 2 ) = 0. 14C 

VXII= 0.350 


LAST STAGE N = 2 WORK. SPLIT* 0.440 EFFICIENCY* 0.901 REACTION* 0.140 
INLET STATIC PRESSURE AT HUB , MEAN, T IP = 679C.849 7788.621 8381.295 

INLET TOTAL PRESSURE* 108 7e.-i77 

INLET TOT VL RELATIVE PRESSURE AT HUB, MEAN, T I P= 8274.052 8718.276 7284. 179 

INLET STATIC TEMPERATURE AT HU T , MEAN, TIP- 2147.631 2216.838 2254.774 

INLET TOTAL TEMPERATURE* 2395.024 

INLET TOTAL RELATIVE TEMPERATURE AT .Hub* MEAN, T IP= 2^48.066 2275.432 2308. 782 

INLET DENSITY AT HUB , KEAN, T I P = C.U59 0.065 0.069 

EXI T STATIC PRESSURE AT HUB ,- M b AN » T I P = 6435.645 6439.276 6441.135 

EXIT TOTAL PRESSURE* 7232.178 

EXIT. TOTAL RELATIVE PRESSURE AT HUB , MEAN* TIP* 8107.205 8585.955 9212.368 

EXIT STATIC TEMPERATURE AT HUB, ME AN, T l P = 2128.581 2128.858 2129.0 1 

EXIT TOTAL T £ MP£R ATURE* 2180.829 


EXIT TOTAL 

RELATIVE 

TEMPERATURE AT HUS, 

MEAN, TIP= 

2245.385 

22 75.389 

2312.764 






EXIT DENSITY AT HUB, 

MEAN, T I P- 

0.056 

0.G66 

. 0.056 








RADE 

DAE 

VAE 

VE VCR 

UEVCR 

WEWCR 

RADI 

BA I 

V A I VIVCR 

IJ I VOR 

« IWOR 

VXI 

VXb 

22. OCO 

-55. 980 

-4.082 

0 ■ 449 o 

0.6323 

0.7794 

21 • 60 9 

-8.972 

52.201 0.6691 

C. 59 35 

0,42 29 

0.410 

>.448 




931.7126 

1310.3090 

1660.8926 



145C..9883 

128 7. .4 . i 

900.412 



18.051 

-51.205 

-4.9 71 

0.4502 

0.5188 

0.7016 

18.051 

18.499 

5 7. T 60 0.7541 

C.4 55 7 

. .4437 

J. »10 

•5, 448 




932.8574 

1075. LO05 

1483.1208 



1635.4931 

1-75.1 ;t4 

*,37.8549 



14.102 

-45.432 

-6.352 

0.4512 

0.4043 

0.6306 

14.493 

43.586 

62.517 0,8886 

C. 3980 

0.5844 

0 . 410 

0.448 




935.0899 

839.9081 

1324. 1620 



1927. 1086 

863,1671 

1227.9766 




FIRST STAGE N= 1 WORK SPLIT* 0.560 EFF IC IENC Y= u,890 RE AC T I ON- 0.180 

INLET STATIC PRESSURE AT HUB , ME AN, T I P = 11119.911 12ol7.624 13629.692 
INLET TOTAL PRESSURE* 18171.970 

INLET TOTAL RELATIVE PRESSURE AT HUB, MEAN, T ! P= 13476.945 14001.377 14633.638 
INLET STATIC TEMPERATURE AT HUB, ME AN, T IP* 2374.282 2444.716 2488.748 

INLET TOTAL TEMPERATURE* 26,0.000 

INLET TOTAL RELATIVE TEMPERATURE AT HUB, ME AN , TIP = 2482.267 2504.280 2530.010 

INLET DENS ITY AT HUB, MEAN, T I P= 0.087 T.,0 96 0.1.01 

EXIT STATIC PRESSURE AT HUB ,MF AN , T I P= 10291.566 10321.834 10339.600 
EXIT TOTAL PRESSURE* 11215.337 

EXIT TOTAL RELATIVE PRESSURE AT HUB » ME AN,T I P= 13029. bCC 13600.248 14305.593 
EXIT STATIC TEMPERATURE AT PUB .MEAN, T IP* 2347.864 2349.460 2350.395 

EXIT TCTAL TEMPERATURE* 2395.024 

EXIT TOTAL RELATIVE TEMPERATURE AT HUB , MEAN, T I P= 2479.576 2504.279 2533.747 

EXIT DENSITY AT HUB ,ME AN, T I P= 0.081 0.081 0.081 


RADE 

BAE 

VAE 

VFVCR 

UE VCR 

WEWCR 

RADI 

BA I 

VAI 

V IVCR 

UIVUR 

WlWuR VXI V XE 

21.226 

-61.336 

-13.531 

0.3774 

0.5829 

0,74 38 

20.849 

12.095 

61,318 

0 . 70 15 

0.54 : 3 

5,35 31 j . 3 3 7 ^.367 




818.5040 

1264.2037 

1659,0319 




1603. 3556 

1241. 7630 

787.01 )2 

18.051 

-58.533 

-15.800 

0.3814 

0.4957 

u.6875 

18,051 

35,534 

64.6 56 

0 • 7866 

0.4 724 

0.4264 j.337 J.3G7 




927.0335 

1075. 1C85 

1524.4901 




1797.7037 

1075. 10 E5 

945.65 *2 

14.876 

-55.532 

-18.951 

0.3880 

0.4086 

u. 6372 

15.253 

52.813 

68,187 

0.9062 

C.397-* 

C. 5767 0. 337 ,.3o7 




341.3916 

866.0133 

1406.1247 




2071.006 2 

908.4540 

1*73.1968 


INITIAL SURFACE AREA CALCULATIONS 


BLADE ENTERING ANGLE, BET A 1 ,FQ ! ROWS 1,2, 

BLADE LEAVING ANGLE, BET A 2, FOR RUhS 1,2.3 
BLADE CAMBER ANGL t , PHI , FOR ROWS 1,2,;>,4= 

BLADE CAMuER LENGTH TiJ CHORD UNGTH RATIO FOR ROWS 
AXIAL CHORDS FOR KQwS 1,2, 3, 4* 2.798 

SINGLE BL^DE AREAS FOR ROWS 1,2. 3, 4= 

PITCH FOR ROWS l ,2 ,3,4* 2.560 

MEAN RADIAL C I RCUMF EREMCF- 113.418 

NUMBER OF BLADES IN ROWS 1,2,3,4- 
TOTAL BLADE SURFACE AREAS IN ROWS 1,2, 3,4= 21C2.041 


(.4 = 


90 

.000 

125 

.534 

105 

.800 

108. 499 

4 = 


25. 

344 

31. 

467 

32. 

940 

38.795 


40, 

.421 

2 7. 

,341 

33. 

469 

24 

.425 

FOR 

ROWS 1, 

2,3.4= 


1 ,49 3 


1.320 

1.328 

t 


L.99 

1 

2.244 

1,501 


47. 

447 


31.859 


40.749 


27.728 


1.546 


1.931 


1.205 




44. 



73. 


59. 


94. 



1.212 
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COOLANT CALCULATION 


INPUT DATA** 


PTD* 0. 

FARC* 

0.023 

NPTD= 12 



T PR =1460. 000 

PRN* 

C * 7 10 

TPR* 1660. 000 

PRN* 

0.7C9 

TPR* I860. COO 

PRN* 

0,708 

TPR=2060 • COO 

PRN* 

0.706 

TPR=2260. >00 

PRN* 

0.705 

TPR*2460. 000 

PRN* 

0.704 

TPR=2660.C00 

PRN* 

0.703 

TPR* 2860. COO 

PRN* 

0.702 

TPR=3060. >00 

PRN* 

0.70C 

TPR=3260. 003 

PRN* 

0.695 

TPR* 3460. >'•00 

PRN* 

0.693 

TPR*3660 . 000 

PRN* 

0.686 

NTS* 1 

TMS=2180.000 

THERME* 0.700 

CPC* 

0.267 


PF= 0.150 

DTY* 400.000 

V I S=0. 240E—D4 

C0N*0. 9506-05 

V I S-0. 260E-04 

CDN=0. 108E-04 

VI S*0 • 280E-04 

CON*0 • 117E-04 

V 1 S=0 . 300E-04 

CON=0 • 127E-04 

V I S=C. 320E-04 

C0N=0 . 137E-04 

V I S=0, 340E-G4 

CON=Q . 147E-04 

V I S*0 • 355E-04 

CON=0. 1606-04 

V I S=0. 370E-04 

C0N=0. 170 E— 04 

V I S = 0 • 390E-C4 

CON*C ■ iaOfc-04 

VI 5*0. 405E-&4 

CON=0. 192E-04 

V I S*0. 420E-04 

C0N*0.205F-04 

VI S=0.435F-04 

C0N=0 . 220 E-04 

NTR= 1 

T MR* 208 C.000 


TNALLA=2060. 000 


INITIAL TOTAL BLAUE SURFACE AREAS FOR ROWS 1,2, 3, 2102*041 2337.120 2393.205 2609.172 

COOLANT INLET TEMPERATURE* 1204.000 

ASSUMED BLADE THERMAL COOLING EFFECTIVENESS* 0.700 

INITIAL TOTAL EXIT TEMPERATURE RELATIVE TO ROWS 1,2, 3,4* 2660.000 2504.279 2395.024 2275.389 

INITIAL TOTAL INL^T TEMPERATURE RELATIVE TO ROWS 1,2, 3, 4= 2660.000 2504.288 2395.024 2275.432 

REVISFD TOTAL INLET TEMPERATURE RELATIVE TO ROWS 1,2,3, 4= 2660.000 2479.326 2354.005 2227.523 

REVISED TOTAL EXlf TEMPERATURE RELATIVE TO ROWS 1,2, 3,4= 2635.038 2463.259 2347.116 2223.450 

REVISED TOTAL RLAuE SURFACE AREAS FOR ROWS 1,2, 3,4= 2050.533 2307.184 2385.141 2606.835 

ASSUMED METAL TEMPERATURE FOR ROW! 11= 2180.000 BLADE COC.LANT FLOW RATIO* 0.016565 

ASSUMED M •„ T A l TEMPERATURE FOR ROW! 21= 2080.000 BLADE COULANT FLOW RATIO* 0.012328 

ASSUMED MCTAL TEMPERATURE FOR ROW! 31= 218C.000 BLADE COOLANT FLOW RATIO* 0.005897 

ASSUMED METAL TEMPERATURE FUR ROW! 41= 2080.000 BLADE COULANT FLOW RATIO* 0.003933 

ASSUMED WALL METAL TEMPERATURE* 2O60.CC0O0O 
WALL CC PLANT FLOW RATIO FDR ROW! 11= 0.015231 

WALL CCDLANT FLOW RATIO FOR ROW! 21= 0, 005888 

WALL CCQLANT FLOW RATIO FOR ROw! 31= 0.CC3C47 

WALL CCOIANT FLOW RATIO FOR R0W1 41= 0.000223 

TURBINE WHISHT FLOW AT EXIT OF ROWS 1,2, 3, 4= 328.818485 334.921719 337.918358 339.310982 

AVERAGE BLADE HEIGHTS FOR ROWS 1.2, 3, 4 WITH BLADE AND WALL COOLANT* 5.546 6.039 6.904 7.726 

REVISED VXI AND VXE FOR FIRST AND SECUND STAGE ROTORS* 0.338 0.370 0.414 0.454 

TURBINE EXIT TOTAL TEMPERATURE WITH BLADE AND WALL CuGLANT* 2112.457 

AVERAGE H.T. COEFFICIENTS FOR ROWS 1,2, 3,4= C.111E+GG 0.106E+00 C.859E-01 0.808E-01 

CRITICAL ABSOLUTE VELOCITIES AT EXITS OF ROWS 1,2, 3, 4= 2275.795 2151.015 2147.865 2048.452 

CRITICAL RELATIVE VELOCITIES AT EXITS OF ROWS 1,2, 3,4= 2207.530 2200,365 2092.429 2090.515 


FAR* C.023 
PT0=18734. CO 
DHT* 1 4 C • 1 J 8 
NR* 2 

WS121* 0.440 

PC* 0.970 
RFI2.J* 0.14. 
VXII* 0. 350 


LAST STAGE N= 2 WORK SPLIT* 0.440 EFFICIENCY* 0.901 REACTION* 0.140 
INLET STAT IC PRESSURE AT HUfi.MEAN, TIP* 6890.863 7817.989 8378.736 

INLET TOTAL PRESSURE* 10953.929 

INLET TOTAL RELATIVE PRESSURE AT HUB , MEAN, TIP* 3391.906 8825.383 9372.423 

INLET STATIC TEMPERATURE AT HUB, MEAN , T IP* 2154,489 2218.339 2254.177 

INLET TOTAL TEMPERATURE* 2398.371 

INLET TOTAL RELATIVE TEMPERATURE AT HUB , ME AN, T I P= 2254.996 2281.426 2313.391 

INLET DENSITY AT HUB .MEAN, T IP* G.059 0.065 0.069 

EXIT STATIC PRESSURE AT HUB , ME AN , T IP* 6523.945 6527.625 6529.510 

EXIT TOTAL PRESSURE* 7329.942 

EXIT TOTAL RELATIVE PRESSURE AT HUB , MEAN, T IP* 8214.524 8698.293 9331.188 

EXIT STATIC TEMPERATURE AT HUB .MEAN, T IP* 2134.497 2134.775 2134.918 

EXIT TOTAL TEMPERATURE* 2192.805 

EXIT TOTAL RELATIVE TEMPERATURE AT HUB, MEAN, T I P= 2251.377 2281.381 2318.756 

EXIT DENSITY AT HUB, MEAN ,T I P= 0.057 0.G57 0.057 


CASE=144 

VELOCITY DIAGRAMS »•»•*•«•, »«••*,«»••• 

INPUT DATA** 

TCI = 1204.G00 CL'RRF* 1.000 CORRFW* UGOO 

Y* 0.0631 YB= 0.0387 TTO*2660.0C0 

RPM=6825 . GOO DRTE* 0.1C0 WO* 321.117 

RTE= 22. COO RHRTE* 0.641 N= 2 

NVUM* I NRTE* 1 WS ! 1 ) * 0.560 

GAMMA* 1.301 CP* G.300 R= 53.980 

ETAT* 0.900 ETA (11= 0.890 REI D * 0.180 

VUMVCR* -0.039 VXEL* 0.800 VXE I = 0.80C 

VXI L= 0.550 


RADE 

BAE 

VAE 

VEVCR 

UEVCR 

WEwCR 

RADI 

BA I 

22.000 

-55.969 

-4.085 

0.4492 

932.1907 

0.6314 

1310.3090 

0.7785 

1661.2326 

21.473 

-8.575 

18.051 

-51 • 193 

-4.975 

0.4498 

933.3380 

0.5181 

1075.1085 

0.7009 

1483.5035 

18.051 

16.666 

14.102 

-45.422 

-6.358 

0.4509 

935.5755 

0.4048 

839.9081 

0.6300 

1324.5949 

14.629 

40.625 


VAI VIVCR Ul VCR WIWCR VXI VXE 

50.684 0.6779 0.5893 0.44.'3 .430 0.448 

1471.2487 1278. 94*1 942.8137 

55.456 0.7575 0.4354 0.4598 0.430 0.448 

1643.9457 1C 75. 1065 973.1559 

60.843 0.8817 0.4J15 0.58370.4300.448 

1913.3829 &7l .26 89 1228.3160 
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FIRST STAGE N= I WORK SPLIT* 0.560 EFFICIENCY* 0.890 REACTION* 0.180 
INLET STATIC PRESSURE AT HUB , MEAN , TI P= 10994.290 12534.069 13571.334 
INLET TOTAL PRESSURE* 18171.979 

INLET TOTAL RELATIVE PRESSURE AT HUB , MEAN , TIP* 13418.820 13947.383 14586.089 
INLET STATIC TEMPERATURE AT HUG. ME AN , T IP* 2368.049 2440.961 2486.279 

INLET TOTAL TEMPERATURE* 2660.000 

INLET TOTAL RELATIVE TEMPERATURE AT HUB, MEAN , T IP* 2479.786 2502.051 2528.105 

INLET DENSITY AT HUB , MEAN, T I P= 0.086 0.095 0.101 

EXIT STATIC PRESSURE AT HUB » ME AN, T IP* 10207.071 10225.395 10236.647 
EXIT TOTAL PRESSURE* 11292.710 

EXIT TOTAL RELATIVE PRESSURE AT HUB, MEAN, T IP* 13032.690 13559.028 14198.323 
EXIT. STATIC TEMPERATURE AT HUB .MEAN, T IP* 2342.935 2343.908 2344.5J4 

EXIT TOTAL TEMPERATURE* 2398.371 


EXIT TOTAL 

RELATIVE 

TEMPERATURE AT HUB, 

MEAN, TIP* 

2479.221 

2502.035 

2528.847 







EXIT DENSITY AT HUB, 

ME AN, TIP* 

0.081 

0.081 

0.081 









RADF 

BAE 

VAE 

VEVCR 

UEVCR 

WEWCR 

RAO I 

BA I 

VAI 

VI VCR 

U I VCR 

WlWCR 

VXI 

VX6 

20.960 

-57.868 

-10.273 

0.4144 

0.5752 

0.7465 

20.883 

12.694 

61.401 

0. 7066 

0.54A2 

0.35 56 

0. 338 

0.408 




899.2326 

1248.3378 

1663.5103 




1614.8746 

1243.7761 

792. 39 . 5 



18.051 

-54.951 

-11.884 

0.4166 

0.4954 

0.6951 

18.051 

36.174 

64.768 

0.7934 

0.4 7 4 

0 . 4320 

0 . 338 

0,408 




904. 1979 

1075.1085 

154U.6920 




1813. 3143 

If. 75. 1085 

957.6276 



15. 142 

-51.787 

-14.083 

0.4203 

0.4156 

0.6483 

15.219 

53.354 

68.331 

0, 9160 

0.3966 

0.5869 

0 . 338 

0.408 




912. 2343 

901.8792 

1430.3347 




2093.4732 

906, 44 .9 

1295.12 51 




INITIAL SURFACE AREA CALCULATIONS ******************** 

BLADE ENTERING ANGLE, BETA 1,F0R ROWS 1,2, 3, 4= 90.000 126.174 101.864 106.666 

BLADE LEAVING ANGLE, BETA 2, FOR ROWS 1,2,3, 4= 25.232 35.049 34.544 38.807 

BLADE CAMhER ANGLE , PH I , FOR ROWS 1,2, 3, 4= 40.297 23.337 34.848 25.207 

BLADE CAMBER LENGTH TO CHORO LENGTH RATIO FOR ROWS 1,2,3, 4* 1.489 1.268 1.329 1.213 

AXIAL CHORDS FOR ROWS 1,2, 3,4= 2.832 1.913 2.110 1.474 

SINGLE BLADE AREAS FOR ROWS 1,2, 3,4* 48.491 28.271 36.055 26.756 

PITCH FOR ROWS l, 2, 3, 4= 2.587 1.336 1.927 1.210 

MEAN RADIAL CIRCUMFERENCE* 113.418 

NUMBER OF BLADES IN ROWS 1,2, 3, 4* 44. 85. 59. 94. 

TOTAL BLADE SURFACE AREAS IN ROWS 1,2, 3, 4= 2125.663 2399.469 2122.615 2507.059 


COOLANT CALCULATION 
INPUT DATA** 


PTD* 0. 


FARC* 

0.023 

PF= C. 150 

DTY* 400.000 

NPTD* 12 






TPR=146Q. }00 


PRN* 

0.710 

VIS*U. 240E-Q A 

CON=0, 950 £-05 

TPR=1660. COO 


PRN* 

0,709 

V I S=0« 260E—04 

C0N=0 • 108 E-04 

TPR=1 860 • COO 


PRN* 

0.708 

V I S= 0 • 2 8 0 E - 0 4 

C0N=0. 117 E-04 

TPR=2060. 000 


PRN* 

0.706 

VIS=0.300E-04 

C0N*0 • 12 7 E-04 

TPR=2260. 000 


PRN* 

0.705 

V I S=0. 320E-04 

CON*0. I37E-0A 

TPR=2A60. r 00 


PRN* 

0.704 

V I S*0. 340E-Q4 

C0N*0. 147E-04 

TPR=2660. 000 


PRN* 

0.703 

V I S=0 .355E-04 

CON*0 . 160 E-04 

TPR=2860. 000 


PRN* 

0.702 

V I S=C. 370E-O4 

C0N=0 . 170 E-04 

TPR* 3060. 300 


PRN* 

0.700 

VI 5=0. 390E-04 

CON=0. 180 E-04 

T PR* 3260 .000 


PRN* 

0.695 

V I S=0. 405E-Q4 

C0N=0. 192 E-04 

TPR=3460.000 


PRN* 

0,693 

VIS*0.420F-04 

C0N*0. 205 E-04 

TPR=3660. 000 


PRN* 

0.686 

VI S=C.435E-04 

C0N=0, 220 E-04 

NTS= 1 


T MS =2180.000 

NTR* 1 

TMR=208 0.000 

THERME* 0.700 

CPC = 

0.267 

THALLA=2060. 000 


INITIAL TOTAL 

BLADE SURFACE 

AREAS 

FOR ROWS 1,2, 3, 4= 

2125.663 2399.469 

2122.615 2507.059 

COOLANT INLET 

TEMPERATURE* 

1204 

.000 




ASSUMED BLADE THERMAL COOLING EFFECTIVENESS- 0.700 

INITIAL TOTAL EXIT TEMPERATURE RELATIVE TO ROWS 1. 2,3*4= 2660.000 2502.035 2398.371 2281.381 

INITIAL TOTAL INLET TEMPERATURE RELATIVE TO ROWS 1,2, 3,4= 2660.000 2502.051 2398.371 2 281.426 

REVISED TOTAL INLET TEMPERATURE RELATIVE TO ROWS 1,2, 3,4= 2660.000 2477.096 2356.921 2233.581 

REVISED TOTAL EXIT TEMPERATURE RELATIVE TO ROWS 1,2, 3,4= 2635.045 2460,591 2350.526 2229.434 

REVISED TUT Air BLADE SURFACE AREAS FOR ROWS 1, 2,3,4= 2068.574 2363.058 2111.156 2498.7A9 

ASSUMED METAL TEMPERATURE FOR ROW l 1 ) = 2180. 000 BLADE COOLANT FLOW RATIO* 0*016708 

ASSUMED METAL TEMPERATURE FOR ROW ( 21 = 2080.000 BLADE COOLANT FLOW RATIO* 0.012804 

ASSUMED METAL TEMPERATURE FOR ROW! 31= 2180.000 BLADE COOLANT FLOW RATIO* 0.005509 

ASSUMED MtTAL TEMPERATURE FOR ROW! 41= 2080.000 BLADE COOLANT FLOW RATIO* 0.004016 

ASSUMED WALL METAL TEMPERATURE* 2060,000000 
WALL CCOLANT FLOW RATIO FOR ROW! 11= 0.015377 

WALL CCOLANT FLOW RATIO FOR ROW! 21= 0.005504 

WALL CCOLANT FLO* RATIO FOR ROWt 31= 0.003049 

WALL CCOLANT FLOW RATIO FOR ROWt 41= 0.000301 

TURBINE WEIGHT FLOW AT EXIT OF ROWS 1,2, 3,4= 331.759995 337.893864 340.761044 342.207428 

AVER AGF BLADE HEIGHTS FOR ROWS 1,2, 3,4 WITH BLADE AND WALL COOLANT* 5.599 5.781 6.478 7.568 

REVISED VXI AND VXF FOR FIRST AND SECOND STAGE ROTORS* 0.340 0.411 0.434 0.453 

TURBINE EXIT TOTAL TEMPERATURE WITH BLADE AND WALL COOLANT* 2118,617 

AVERAGE H.T. COEFFICIENTS FCR ROWS 1,2, 3,4= 0.1 HE +00 0.108E+00 J.897E-01 0.827E-01 

CRITICAL ABSOLUTP VELOCITIES AT EXITS OF ROWS 1,2, 3,4= 2275.798 2152.347 2149.425 2051.314 

CRITICAL RELATIVE VELOCITIES AT EXITS OF ROWS 1,2, 3,4= 2206.536 2199.173 2095.273 2093.327 


THE LAST SET OF OUTPUT CONTAINS THE CONVERGED SOLUTION 
*01* EXIT IN COOL 
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APPENDIX B 


PROGRAM LISTING 


COMMON /XLOOP/ RY,RYB,RWO,ROHT,NWC, ITT, TESTY 

COMMON /CFCT / TC I ,CORRFW ,CORRF , YW, TPPMTP ( 100 > ,OHC , WPSTAT ( 5 ) ,HTEMP ( 
15) 

COMMON /SV/ XVAI ( 5 ) • XBAI (5) ,XBAE(5) , XVAE ( 5 ) , SA VE » BHS ( 5 ) , 6HR ( 5 ) , RAE 
1(20) » TGA ( 20 ) »WF,RAT(3) » R AH{ 3 ) » RHOI ( 10 ) , RHOE ( 10 ) , VEU 10 ) , VI L 1 10 ) , PP 
2(10) ,TT( 10 ) » CONS T2» CONST l» VI LC I(10)»GAMMA,R»DTR1»DTR2»FAR,WETO»VXI 
32»VXE2»VXI1,VXE1 
COMMON VXIL 

DIMENSION VUK20), VUEI20), VI(20), VEI20), UI(20), UE(20), WUI(20 
1), WUE ( 20 ) » W 1 2 ( 20 ) » WE2 (20) » TRTE ( 20 ) » WCRE(20), BAE(20>, BA I( 20 ) 
2, VAE120), V A I ( 20 ) » VEVCRI20), VIVCR(20), UEVCR120), UIVCRI20), WE 
3WCR ( 20 ) « WIWCR120), TRTH20), WCRM20), DH(20), WS(20), ETA120), R 
4AI(20), RE ( 20 ) * VUMVCR(20>, TETI20), TPE(20), RT ( 50 ) , SVEI20), SWE 
52(20), SV I ( 20) , SW 12(20) » DUMMY(IOO), NDUMMY(iuO), PRTK100), PRTE 
6(100), TSTATE ( 100 ) , TSTATI ( 100) , PSTATE(IOO), PSTATI(IOO), RHOSE ( 1 
700), RHOSI(IOO) 

DIMENSION H I R ( 5 ) , HE R ( 5 ) , CVAE(IOO), TPPCU(IOO), TPPCV(IOO), TPPCV 
1U ( 50 ) , DUM ( 10) tNDUM(lO) 

1 READ (5,51) NCASE 

READ (5,52) TCI , CORRF ,CORRFW 
READ (5,53) FAR , Y , YB 

READ(5,54)TT0,PT0,RPM,DRTE,W0,DHT,RTE,RHRTE,NWC 

PT0=PT0/144. 

N=2 

NR=2 

NVUM=1 

NRTE=1 

REAO (5,55) (WS(I),I=1,N) 

READ (5,56) GAMMA, R, PC 
CP=(GAMMA*R) /(GAMMA-1. ) 

CP=CP/778. 

READ (5,55) ETAT, ( ETA (I) , I=1,N ) 

READ (5,55) (RE ( I ),I*1,N) 

READ (5,55) ( VUMVCR( I ) , I =1 , N VUM ) 

READ (5,55) VXEL, VXE I , VX I I , VXIL 

IF(NWC.EQ.l) Y=YB 

DUM ( l ) =RTE 

OUM (2) =RHRTE 

NDUM ( 1 ) =N 

DUM ( 3) =PT0 

ITT-0 

2221 CONTINUE 
ITT=ITT+1 
WRITEI6, 2223) ITT 

2223 FORMAT ( 1H , 1 6H ITERAT I ON NUMBER, 14) 

IF C ITT . EO . 1 ) GO TO 5000 

Y=R Y 

YB=RYB 

IF(NWC.EQ.I) Y=YB 

WO=RWO 

DHT=ROHT 
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o o 


RTE = OUM(l ) 

RHR TE=DUM { 2 ) 

N=NDUM ( 1 > 

PTO-DUMt 3 ) 

5000 CONTINUE 

WRITE (6,57) NCASE 
WRITE (6,40) 

WRITE (6,50) 

WRITE (6,41) TCI »CORRF,CORRFW»FAR 
TESTY* V 
PT0=PT0*144 . 

WRITE (6,42) Y,YB,TTO,PTO 
PT0=PT0/144» 

WRITE (6,43) R PM , DRT E , WO , DHT 
WRITE (6,44) RTE,RHRTE,N,NR 
WRITE (6,45) NVUM, NRTE,WS( l ) , WS( 2 ) 

WRITE (6,46) GAMMA, CP ,R, PC 
WRITE (6,47) ETAT ,ET A ( l ) ,RE( 1) , RE ( 2 ) 

WRITE (6,48) VUM VCR ( 1 ) ,VXEL, VXEI ,VXI I 
WRITE (6,49) VXIL 

0HC8=l.+.5*( YB/( ( l.-Y) *( l.+FAR ) ) )*(TCI/TTO) 
DHC=DHT*< ( l.-Y ) * ( l. + FAR) )*DHC8 
WETO=WO/( l ,+FAR* ( 1 Y )-Y ) 

WCBTO* YB+WETO 
WF = W0+WCBTCJ/2. 

YW=Y-YB 

MINE=1 FOR SECOND STAGE 
MINE=2 FOR FIRST STAGE 
MI NE=0 

DUMMY! 1 )=TTO 
DUMMY( 2 ) =PTO 
DUMMY( 3 ) *RPM 
DUMMY ( 4) -DRTE 
DUMMY (5 )*WF 
DUMMY(6)*DHT 
DUMMY! 7 ) *RTE 
DUMMY ( 8 ) * RHRTE 
NDUMMY ( 9 ) = N 
NDUMMY ( 10 ) =NR 
NDUMMY til) =NVUM 
NDUMMY ( 12) =NRTE 
DUMMY( 1 3 ) *WS ( 1 ) 

OUMMY( 14) =WS (2) 

DUMMY! 15 ) -GAMMA 
DUMMY ( 16 ) =CP 
DUMMY ( 1 7 ) *R 
DUMMY( 18)=PC 
DUM MY ( 19) = E T A ( 1 ) 

DUMMY ( 20 ) =ET A ( 2 ) 

DUMMY! 21 >*RE(1) 

DUMMY! 22 ) -RE (2 ) , 

DUMMY! 23) - VUM VCR ( 1 ) 

.DUMMY (24)=VXEL 
DUMMY! 25 ) =VXEI 
DUMMY ( 26 ) =VX 1 1 
C J= 778 . 

G = 32. 17 
2222 CONTINUE 
NA=N 
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T TO = DUMMY ( 1 ) 

PT0=0UMMY(2) 

R PM-DUMMY { 3 ) 

DRTE=DUMMY<4) 

WF=DUMMY( 5) 

DHT = DUMMY (6 ) 

RTE -DUMMY ( 7 ) 

RHRTE=DUMMY(8) 
fs = NDUMMY(9) 

NR = NDUMMY( 10 > 

NVUM=NOUMMY( 11) 

NRTE=NDUMMY( 12) 

WS { l )=DUMMY( 13) 

WS(2)=DUMMY( 14) 

GAMMA=DUMMY( 15) 

CP=DUMMY( 16) 

R— DUMMY (17) 

PC=DUMMY( 18) 

ETA( l) = DUMMY! 19) 

ETA(2)=DUMMY(20) 

RE ( l )=DUMMY { 21 ) 

RE ( 2 )- DUMMY ( 22 ) 

VUMVCR( 1)=DUMMY(23) 

VXEL=DUMMY ( 24) 

VXEI=DUMMY (25) 

VXI I=DUMMY{26) 

RT( i)=RTE 
DO 2 KL=2,NRTE 
RT(KL)=RT(KL-1 )+DRTE 

2 CONTINUE 
DHT=DHT*C J 
PTO=PTO* 144. 0 

IF (N.EQ.l) GO TO 4 

PTE=PTO* ( 1 .O-OHT/ ( CP*TTG*ETAT *CJ ) ) ** (GAMMA/ ( GAMMA- 1 .0 ) ) 

PT=PTO 

TT=TTO 

NS=N-1 

DO 3 11=1, NS 
DH( 1 1 ) =WS ( II )*DHT 

TPE( II )=PT*( l.O-DHU I )/(CP*TT*ETA( II )*CJ) ) **( GAMMA/ ( GAMMA-1 .0 ) ) 
PT = TPE( 1 1 ) 

TET t 1 1 )=TT-DH( II )/(CP*CJ ) 

3 TT=TET (II) 

ETA(N)=WS(N) *DHT/ (GP*TT* ( 1 .0— ( PTE/PT ) **.( (GAMMA-1.0 ) /GAMMA) )*CJ) 

4 00 39 KL* 1 , NRTE 
RTE=RT ( KL ) 

RTE=RTE/12.0 

DO 38 K=1 , NVUM 
VUEVCR* VUMVCR (K ) 

KODE=D 

N=NA 

KRTE=0 

C LAST STAGE EXIT 

DH(N)=WS(N)*DHT 
TTE=TTO-DHT / ( CP*C J ) 

PTE—PTO* ( 1 . O-DHT/ ( CP* TTO*ETAT*C J ) )**( GAMMA/ ( GAMMA- 1.0) ) 

TTE2=TTE 

PTE2=PTE 

RAPS=3.1416#RPM/30.0 
KCC = 1 
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WHIC=0.0 
LTE =RTE*R APS 
RHE=RTE*RHRTE 
RME=(RHE+RTE>/2.0 

VUEf‘l=V l JE VCR* SORT l 2 »0*G*G AMM A*R*TTE/ ( GAMMA* 1 . Q ) } 

CE=VUEM*RME 
VUE ( 1 ) =CE/RT E 
VUE ( 2) =CE/RME 
VUE ( 3) =CE/RHE 

AE=3.1416*RTE**2*(1.0-RHRTE**2) 

CALL DENR ( VXE .G.GAMM A,CP, R, TTE , PTE , VUE ( 2) , WF , AE , RHOM ) 

IF ( VXE.GT.VXEL*SQRT { 2. 0 *GAMMA*G*R*TTE/ ( GAMMA+1 .0 ) >1 GO TO 36 
CALL VXR (VUE, VE, VXE, WF, TTE, PTE, G, GAMMA, CP, R,RTE»RME,RHE) 

UE ( 1 ) =UTE 
UE(2 )=RME*RAPS 
UE ( 3 ) = RHE*RAPS 
ViUE ( 3) =VUE ( 3 ) -UE 1 3 ) 

WE2(3)=VXE**2+WUE(3)**2 
LAST STAGE INLET 
mWEH2=l.-RE(N) 

CI=(G*DH(N)+RAPS*CE) /RAPS 
TTI=TTE+DH(N)/(CP»CJ ) 

PTI =PTE/ ( ( l.O-DH(N)/ (CP*CJ*TTI*ETA(N) ) )**( GAMMA/ ( GAMMA-1 .0 ) ) )*PC 
TT I 2=TT I 
PTI2=PTI 
RM I =RME 

RHI2=( RME-RHE ) *Q.55+RHE 

RHRT I 1=RHRTE- .02 

RHRT I2=l.0/(2.0*RMI/RHI2— 1.01 

CALL RHRT2 ( RHRT 1 2 , GAMMA ,CP , G»R , T T I , PT I , WF , RM I ,C I , *2222 ) 

RHRTIL—RHRTI 2 

KCHK=0 

TOLER=. 0002 

I ND = 1 

AI = 3. 14 16*4. 3*RM I **2 * ( 1. O-RHRT 1 1 )/ ( 1 .0+RHRTI 1 ) 

RTI=SQRT(AI/ (3.1416* ( 1.0-RHRTI 1**2 ) ) ) 

RHI =RT I *RHRT 1 1 
VUI (1)=CI/RTI 
VUI (2)=CT/RMI 
VUI ( 3) “Cl /RH I 

CALL DFNR ( V XI ,G , G AMM A ,C P, R , TTI , PT I , VU I ( 2 ) , WF , AI , RHOM > 

CALL VXR ( VU I ,VI ,VXI ,WF, TTI , PTI, G, GAMMA, CP, R,RTI,RMI, RHI ) 

U I ( 3 1 =RHI *RAPS 

WH I 2=( VUI (3) -U 1(31 )**2+VXI**2 
WHI 2=-WHI 2 

CALL MAXIM ( RHRT 1 1 ,RHRT I 2, WH 12 , IND , TOLER ) 

IF ( INO.LT, 61 GO TO 6 
WHl 2-- WHI 2 
RHRTI2=RHRTE— .002 
WC=WIWEH2*WE2(3) 

A 11=3. 1416*4. 0*RM 1**2* (1.0-RHRT 12)/ ( U0+RHRTI2> 

RTI=SQRT( All/ (3.1416*1 1. -RHRT 12**21) ) 

RHI =RTI *RHRT 12 
VUI ( 1)=CI/RTI 
VUI (2)=CI/RMI 
VUI ( 3 ) =C I /RH I 

CALL DENR ( V XI ,G , GAMMA ,CP, R , TT I , PT I , VU I ( 2 ) , WF , A 1 1 , RHOM ) 

CALL VXR (VUI.VI.VXI ,WF, TTI, PTI, G, GAMMA, CP, R,RTI,RMI, RHI ) 

UI ( 3)=RHI*RAPS 



WHI21=( VUI (3J-UI <3>)**2+VXI**2 
IF (WHI2l.GE.WC) GO TO 8 
IF (WHI21.LE.WC) GO TO 9 

8 IF (WHI2.LE.WC) GO TO 12 
IF (KCHK.EG.l) GO TO 10 
RHRT I 2*RHRT I L 
WHI2E=WHI21 

KCHK=1 
GO TO 7 

9 IF ( WHI 2.GE.WC ) GO TO 12 

IF (KCHK.EO.l) GO TO 10 
RHR TI2=RHRTIL 
WHI2E=WHI 21 

KCHK^l 
GO TO 7 

10 WH 12 1=WHI 2E 

IF (KCC.EQ.l) WCC=WC 

IF (AB.S(WHI21-WC) .GT.ABSt WHI C-WCC) ) GO TO 11 

RTE-RTE-. 125/12.0 

WH I C=WB 121 

WCC^WC 

KCC^O 

GO TO 5 

11 RTE=RTE+ 1.0/ 12.0 
WH I C=0 • 0 

KCC-1 

wcc=wc 

KRTt-KRTE+l 

IF (KRTE.GT.30) GO TO 36 
GO TO 5 

12 CALL ROOT ( WHI 2 « WHI 2 1 » RHRT I I *RHRTI 2 1 WC » RHRTI » TT I , PTI f GAMMA * CP* WF 9 R 
1 f VUI » V! , RT I *RM 1 1 RHI * RAPS *C t • VX1 ) 

U! m*RTI*RAPS 
UI ( 2 ) =RMI *R APS 
U 1 1 3 ) = RH I *RAPS 
AR- ( RT I /RTE ) **2 
GO TO 28 

C FIRST STAGE 

13 DH { N ) = W S t N ) * OH T 
WIWEH2-1.-RE ( N ) 

RT I R-RT I 
ARLS=AR 
PTE=PTI/PC 
TTE=TT I 

TTI=TTE+DH(N)/(CP*CJ ) 

PTI=PTE/( t l.-D.HiNI/tCP*TtI*CJ*ETA(N) I )**■{ GAMMA/ ( GAMMA-1 . 0 ) ) ) *PC 

PTE l-PTE 

TTE i=TTE 

T T 1 1 = T T I 

PTII^PTT 

C FIRST STAGE EXIT 

RME =RM I 

C FIRST ASSUMED MEAN EXIT WHIRL 

VUEM=0. 1 * SORT l 2 »0*G*G AMM A*R*TTE/ ( GAMMA* 1.0 ) ) 

VUEa^VUEM 
GO TO 15 

C SECOND ASSUMED MEAN EXIT WHIRL 

LA VU£M= ( - # A ) *SQRT ( 2 • 0*G*GAMM A*R*TTE/ ( GAMMA + 1 • 0 ) ) 
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15 


CE- VUEM*RME 
AR=ARLS 
GO TO 17 

16 IF (AR.EQ.1.1) GO TO 37 
AR= AR+ . 01 

IF (AR.GT.1.1) AR= 1.1 

17 RTE=RTIR*SORT(AR) 

RHE=2.0*RME-RTE 
AE=3.1A16*(RTE**2-RHE**2 I 
VUE ( 1 ) =CE/RT E 

VUE (2) =VUEM 
VUE ( 3) =CE/RHE 
UE ( 1 ) =RTE*RAP$ 

UE( 2)=RME*RAPS 
UE(3)=RHE*RAP$ 

IF IABS ( VUE 1 1 ) ) .LT.0.5) VUE (11 *0.0 
IF (ABS(VUE(2) ). LT.0.5) VU6(2)=0.G 
IF (ABS(VUE(3) ). LT.0.5) VUE(3)=0.0 

CALL DENR ( VXE ,G .GAMMA ,CP, R , TTE, PTE , VUE l 2 ) , WF , AE.RHOM ) 

IF (VXE.GT.VXEI*SQRT(2.0*GAMMA*G*R*TTE/(GAMMA+1.0> ) » GO TO 16 
CALL VXR ( VUE t VE * VXE » WF. TTE » PTE » G» GAMMA » CP » R ,R TE » RMF. » RHE ) 
WUE(3)=VUE(3)-UE(3) 

WE2(3)=WUE(3)**2+VXE**2 
C FIRST STAGE INLET 

CI=(G*DH( N)+RAPS*CE) /RAPS 
RMI =RME 
GO TO 19 

18 IF (AR.EQ.1.1) GO TO 37 

AR= AR+. 01 

IF (AR.GT.1.1) AR=1. I 

19 RTI =RTE*SQRT ( AR ) 

RHI=2.D*RMl-RTI 

AI = 3.1A16# ( R T I **2-RH 1**2 ) 

VUI ( 1J=CI/RTI 
VUI (2)=CI/RMI 
VUI (3)=CI/RHI 

CALL DENR ( V X I , G, GAMM A , C P, R , TT I , PT l , VUI ( 2 ) . WF , A I , RHOM ) 

IF ( VXI » GT . V X I I »SQRT (2.0 *GAMMA*G*R*TTI / ( GAMM A + 1.0 ) ) ) GO TO 18 
CALL VXR ( VUI ,VI ,VXI .WF.TTI, PTI , G , GAMMA, CP , R , R T I .RMI.RHl ) 

UI ( 1 )=RTI *RAPS 
UI (2 )=RMI*RAPS 
U I ( 3 ) =RHI *RAPS 
WUI (31 =VU I ( 3 )— UI (3) 

WI2(3)=WUI ( 3 ) **2+VXl **2 
WlWb2l = Wl2(3J,/WE2(3) 

IF (KODE.EQ. 1) GO TO 26 
IF (VUEM.EQ. VUEA) GO TO 20 
GO TO 21 

20 WlWE2A=WIWE2l 
VUEMA=VUEM 

GO TO 1A . 

21 WIWE22=WIWE21 



IF 

(WIWE2A.GT.WIWEH2 ) 

GO 

TO 

22 


IF 

(WIWE2A.EQ.WIWEH2 ) 

GO 

TO 

28 


IF 

(WIW62A.LT.WIWEH2 ) 

GO 

TO 

23 

22 

IF 

(WIWE22.LT. W1WEH2 ) 

GO 

TO 

24 


IF 

(WIWE22.EQ.WIWEH2 ) 

GO 

TO 

28 


GO 

TO 35 




23 

IF 

(WIWE22.GT.WIWEH2 ) 

GO 

TO 

24 
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IF (WIWE22.EQ.WIWEH2) GO TO 28 
GO TO 35 

24 HA= VUEM-^VUEMA 
HAN=HA 
VUEA=VUEMA 

25 AB=WIWEH2-WIWE2A 
CB=WlWE22-WlWfc2A 

IF (ABS (HAN) .LT. 0.0001) HAN=0.0 

H=AB/CB*HAN 

VUEM1=VUEA+H 

VUEM=VUEM 1 

KOD£=l 

GO TO 15 

26 IF ( ABS(WIWE21-WIWEH2).LE. ( ABS ( WI WE22-W l WE2A) *.000 i > ) GO TO 28 

IF (AB*(WIWE21-WIWEH2).LT.0.0) GO TO 27 

HAN = H 

V» I WE22 = WI WE2 1 
GO TO 25 

27 HAN=HAN-H 

V* I WE2A=W I WE2 1 
VUE A=VUEM 1 
GO TO 25 

28 MINE=MI NE+1 

VCRE=SQRT(2.*G*GAMMA*R*TTE/(GAMMA+i. ) ) 

VCR I =S 3RT { 2 . 0*G*GAMM A *R* TT l / ( GAMMA* 1 .0 ) ) 

VXI V=VX I / VCR I 
VXEV=VXE/VCRE 
RAE ( 1) =RTE*12.0 
RAE(2)=RME*12.0 
RAE ( 3 ) =RHE*12 .0 
RAI( 1) =RTI*12.0 
RAI (2)=RM 1*12.0 
RA I ( 3 ) =RH 1*12.0 

IF (MINE.EQ.l) HER (2 ) =RAE( 1 ) -RAE I 3 ) 

IF I MINE.EQ.l) HIRI2 ) =RA It 1 ) -RAI ( 3 ) 

IF (MINE.EQ.2) HER (1) =R A El 1 ) -R AE ( 3 ) 

IF (MINE. EQ. 2 ) H IR (1) =RA I ( 1 ) -RA I ( 3 ) 

IF (MINE.EQ.2) R AT ( 1 ) =RA I C 1 » 

IF (MINE.EQ.2) RAH ( 1 ) =RA I ( 3 ) 

BHS ( 1) =HI R ( 1 ) 

BhfR ( 1 ) = ( HER ( 1 ) +H IR ( 1 ) )/2. 

BHS ( 2) = ( H IR ( 2 ) +HER ( 1 ) )/2. 

BHR (2 ) = ( HER ( 2 ) +HIR (2 ) ) /2 . 

DO 29 1=1,3 

WUE ( I ) = VUE ( I )-UE(I ) 

WE2( I) =WUE( I )**2+VXE**2 
E AE ( I ) = AT AN2 ( WUE ( I ) » VXE ) *57. 296 
VAE( I ) =ATAN2 ( VUE ( I ) , VXE ) *57.296 
VEVCR( I )=SQRT ( VE( I ) ) / VCRE 
UEVCR( I )=UE( I ) /VCRE 

TRTE ( I )=TTE-VE( I ) / ( 2 . 0*G*C j*CP ) +WE2 (I) / ( 2.0*G*CJ*CP) 

WCRE ( I ) =SQRT ( 2.0*G*G AMMA*R*TRTE ( I ) / ( GAMMA+1.0 ) ) 

WEWCRU ) = S0RT(WE2( I) ) /WCRE ( I ) 

wun i ) =vui ( i ) -u i ( i ) 

WI2 ( I ) =WU I l I )**2+VXI **2 
BA I ( I ) = AT AN2 { WUI 1 1 ) , VX I ) *57 . 296 
V A I l I ) = AT AN2 ( VUI ( I ), VXI ) *57. 296 
V I VCR ( I )=SQRT (VI ( I ) ) /VCR I 
U I VC R ( I )=UI ( I ) /VCRI 
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TRTI (I )=TTI-VI ( I )/(2.0*G*CJ*CP)+WI2( I ) / ( 2 . 0*G*CJ*CP ) 

PRTI (l )=PTI*( ( TRTI ( I ) /TT I ) ** (GAMMA/ (GAMMA- 1 . ) ) ) 

PRTE ( I ) = PTE* I ( TRTE ( I ) /TTE) ** (GAMMA/ (GAMMA-1 .) ) ) 

TSTATE ( I) =TTE* ( 1 .-( ( GAMMA-1. )/(GAMMA+l. > )*VEVCR( I )**2. 1 
TSTATI ( I) =TT I* (l.-{( GAMMA-1. )/(GAMMA+l. ) ) *V I VCR ( I ) **2 . ) 

PST AT I ( n-Pri* (!.-(( GAMMA- 1. )/(GAMMA+l. ) )*VIVCR( I » **2. ) ** ( GAMMA/ ( G 
1AMMA-1. ) ) 

PST ATE ( I )=PTE*(l.-( (GAMMA-1. > / ( GAMMA+1 . ) ) # VEVCRt 11**2. )** (GAMMA/ (G 
1AMMA-1. ) ) 

RHOSE( I )=PSTATE( I )/(R*TSTATE( II) 

RHOSI( I )=PSTAT I ! I )/( R*TSTATI (I)) 

WCRI (I )=SQRT(2.0*G*GAMMA*R*TRTI (I )/( GAMMA+1 . 0) ) 

WIWCRl I i - SORT l W 1 2 ( I ) ) /WCRI II ) 

SVE ( I ) -SORT ( VE ( I)') 

SWE2 ( I ) =SQRT ( WE2 ( I ) ) 

SVI ( I ) * SO R T ( V I ( I ) ) 

29 SWI 2 ( I ) *SQRT ( W.I.2 ( I ) ) 

XBAE ( M I NE ) *8AE ( 2 ) 

XVAE (MINE ) =VAE ( 2 ) 

XBAI (MINE>=BAI(2) 

XV A I (MI N E ) = V A I ( 2 ) 

IF (MIME .EG. 2) TGA(l)=TTI 
IF (MINE.EQ.2) TGA(2I=TRTI (2) 

IF (MINE.EQ.2) CONST 2= < W 12 ( 2 >-VI ( 2 ) ) / ( 2 . *G*C J*CP ) 

IF (MIME.EQ.2) TGA(3 ) =TTE 
IF (MINE.EQ. 1) TGA(4 )=TRTI (2) 

IF (MINE.EQ. 1) CONST 1* (W 12 (2 )-VI ( 21 ) / ( 2. *G*C J*CP ) 

IF (N.FQ.NA) GO TO 30 
GO TO 31 

30 CONTINUE 

WRITE (6,58) N , WS ( N) » ETA <N ) * RE ( N ) 

WRITE (6,63) PSTAT I ( 3 ) » P ST AT I ( 2 ) , PST AT 1 1 1 ) 

WPSTAT (3)*PSTATI (2) 

WRITE (6,71) PTI 
PP(4)=PTI 

WRITE (6,66) PRT I ( 3) , PRT II 2 ) »PRT I ( 1 ) 

WRITE (6,64) TSTAT H 3 ) , TSTAT I ( 2 ) , TSTATH 1 ) 

WRITE (6,70) TTI 
HTEMPI 3)=TT1 
TT( 4 ) =TTI 

WRITE (6,68) TRT I ( 3) » TRT I ( 2 ) , TRTI ( 1 ) 

RHOI ( 4 ) =RHOS 1(2) 

RHOE (3.) =RHOS I ( 2 ) 

RHOE ( 4 ) =RHOSE (2 ) 

WRITE (6,60) RHOS I l 3 ) »RHOS I ( 2 ) « RHOSI { 1 ) 

WRITE (6,62) PST ATE( 3 ) , P ST ATE ( 2 ) .PSTATE ( 1 ) 

WPSTAT(4)*PSTATE(2) 

WRITE (6,73) PTE 
PP(5)=PTE 

WRITE (6,67) PRTE(3),PRTE(2),PRTE(1) 

WRITE (6,65) TSTATE( 3 ) ,TSTATE(2) , TSTATE'C 1) 

WRITE (6,72) TTE 
DTR2=TTE-TT I 
TT ( 5 )=TTE 

WRITE (6,69) TRTE ( 3) , TRTE(2 ) , TRTE ( 1 ) 

HTEMP(4)=TRTE(2) 

WRITE (6,61) RHOSE ( 3 > , RHGSE ( 2 ) , RHOSE ( 1 ) 

GO TO 32 
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31 WRITE (6,59) NNfWS(NN) , E TA (NN ) , RE ( N) 

WRITE (6,6?) PSTATI(3)»PSTATI(2)»PSTATI(1) 

WPS TAT ( 1 ) = PS TAT I (2) 

WRITE (6,71) PTI 
P P ( 2 ) = P T I 

WRITE (6,66) PRTI (3) ,PRTI(2)»PRTI ( 1) 

WRITE (6,64) TSTATI ( 3 ) » TST AT I( 2 ) » TST ATI ( 1 ) 

TT ( 2 ) = T T I 
WRITE (6,70) TTI 
HTEMP( 1 )=TTI 

WRITE (6,68) TRT 1(3) , TRT I < 2 ) , TRT I ( 1 ) 

RHOE ( 1 ) =RHOS 1(2) 

RHOI (2) -RHOS I (2) 

RHO I ( 3 ) =RHOSE ( 2 ) 

RHOE ( 2 ) = RHOS E (2 ) 

WRITE (6,60) RHOS I <31* RHOS 1(2) »RHQS 1(1) 

WRITE (6,62) P STATE! 3 ) , P STATE ( 2 ) , PST ATE ( 1 ) 

WPSTAT (2)=PSTATEl 2) 

WRITE (6,73) PTE 
PP( 3 )=PTE 

WRITE (6,67) PRTEI3) ,PRTE(2) ,PRTE( 1) 

WRITE (6,65) TSTATEI 3 ) , T STATE ( 2.) » T STATE ( 1 ) 

WRITE (6,72) TTE 
DTRi=TTE-TTI 
T T ( 3 ) = T T E 

WRITE (6,69) TRTE (3) , TRTE( 2 ) , TRTE ( 1 ) 

HTEMPl 2 ) =TRTE ( 2 ) 

WRITE (6,61) RHOSE ( 3 ) , RHOS E ( 2 ) , RHOSE ( 1 ) 

32 WRITE (6,74) 

IF (MIME.EQ.I) VX I 2= VX I V *VCR I 
IF (MINE.EO.l) V XE 2= V X E V * VC R E 
IF (MIME.EQ.2) VX 1 1= VX IV*VCR I 
IF (MINE.E0.2) VXE1=VXEV*VCRE 
DO 33 1=1,3 

WRITE (6,75) RAEl I),BAE( I) , V AE ( I ) , VE VCR ( I ) ,UEVCR( I ) , WEWCR ( I ) ,RAI ( I 
1 > ,BAI( I ) , VAI (I) , VI VCR ( I ) , U I VCR ( I ) , WI WCR 1 1 ) , VXI V, VXEV 
WRITE (6,76) SVE ( I ) , U E ( I ) , SWE2 ! I ) , SVI ( I ) «U 1 ( I ) , SWI 2( I ) 

33 CONTIMIJE 

IF (MIME. 60. 21 VEL ( I ) =SV K 2 ) 

IF (MIME.EQ.2) V I L ( 2 ) =SW 12 ( 2 ) 

IF (MINE. EO. 2) VEL ( 2 ) =SW E2 ( 2 ) 

IF (MINE.EQ.2) V I L ( 3 ) =SV E ( 2 ) 

IF (MINE.EO.l) VEL (3 ) =SV 1(2) 

IF (MINE.EO.l) V IL (4 ) =SW 12 ( 2 ) 

IF (MINE.EO.l) VEL (4 ) =SWE2 (2 ) 

IF (MINE.EQ.2) V I LCI ( 2 ) = SV I ( 2 ) 

IF (MINE.EQ.2) VI LCI ( 3 )=SVE( 2) 

IF (MINE.EO.l) V ILC I ( 4 ) =SV I ( 2 ) 

IF (MINE.EO.l) V ILC I ( 5 ) = SVE ( 2 ) 

IF (MINE.EO.l) CVAE( 5 ) =VAE ( 2 ) *3, 1416/180, 

IF (MINE.EQ.2) CVAE( 3 ) =VAE { 2 )*3. 1416/180. 

IF (MINE.E0.2) TPPCU ( 3 ) =UE ( 2 ) 

IF (MINE.EO.l) TPPCU ( 5 ) =UE ( 2 ) 

TP PC V ( 3 ) = V I LC I ( 3 ) 

TPPCV(5)=VILCI (5) 

GRAT=- ( GAMMA-1 . ) / ( 2. *GAMMA*G*R ) 

TPPCVlH 3 ) =TPPC V ( 3 ) *S I N ( C VAE ( 3 ) ) 

TPPCVU ( 5 ) = TPPCV ( 5 ) *S I N (C VAE ( 5 ) ) 

TPPMTP ( 3 ) =GR AT* ( 2 . *TPPCU ( 3 ) *TPPCVU ( 3 ) -TPPCU ( 3 ) *TPPCU( 3 ) ) 
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TPPMTPf 5)=GRAT*( 2.*TPPCU 15 )*TPPCVU(5 )-TPPCU( 5) *TPPCU(5 ) ) 

AN-N-1 

IF (NN.GT.O) GO TO 34 
GO TO 38 

34 N=N-1 

GO TO 13 

35 WRITE (6,77) 

GO TO 38 

36 WRITE (6,78) 

GO TO 38 

37 WRITE (6,79) 

38 CONTINUE 

39 CONTINUE 
CALL SURC 
GO TO 2221 

C 

C 

40 FORMAT (1H ,38HVELQCITY DIAGRAMS ********************,/) 

41 FORMAT (1H , 4HTC I * ,F 8 . 3, 17X , 6HC0RRF= , F8 . 3 , 1 5X , 7HC0RRFW=, F8 . 3 , 1 4X , 4 
1HFAR=, F 8 . 3 ) 

42 FORMAT ( 1H ,2HY=,F8.4, 19 X, 3H YB= , F8 .4 , 18X , 4HTT0=, F8 . 3 , l 7X , 4HPT0= , F8 
1.21 

43 FORMAT (1H , 4HRPM= ,F 8 . 3, 17X , 5HDRTE= , F8 . 3 , 16 X , 3HW0= , F8. 3 , 1 8X , 4HDHT= 

1 » F8 . 3 ) 

44 FORMAT (1H , 4HRTE= ,F 8 . 3 , 17X , 6HRHRTE = , F8 . 3 , 1 5X , 2HN= , I 3 , 24X , 3HNR = , I 3 
1 ) 

45 FORMAT 11H , 5HNVUM= , I 3 , 2 IX , 5HNRTE= , 13, 2 IX , 6HWS (1) = , F 8. 3 , 15X , 6HWS ( 2 
1 ) = » F8. 3 ) 

46 FORMAT (1H , 6HGAMMA® , F8. 3, 1 5X» 3HCP= »F8 . 3 » 18X , 2HR* » F8. 3 , 1 9X» 3HPC= » F 

18.3) 

47 FORMAT (1H , 5HET AT=, F8 . 3 , 1 6X , 7HETA ( 1 ) = , F 8 . 3 , 14X , 6HRE (1) = . F 8 . 3 , I 5X , 
16HRE (2)=,F8.3) 

48 FORMAT (1H , 7HVUMVCR * , F8 . 3 , 1 4X , 5HVXEL=, F8 . 3 , 16X , 5HVXE I = , F 8. 3 , 16X , 5 
IHVXII=,F8.3) 

49 FORMAT (1H , 5HVX I L = , F 8. 3 ) 

50 FORMAT (1H , 12HINPUT DATA**) 

51 FORMAT (13) 

52 FORMAT (F12.6.2F6.4) 

53 FORMAT (3F12.6) 

54 FORMAT (8F6. 1,412) 

55 FORMAT (12F6.1) 

56 FORMAT (4F6.3) 

57 FORMAT ( 1H1, 20X, 5HCASE=» 13, / ) 

58 FORMAT (1H0,14H LAST STAGE N=, 12, 3X, 11HW0RK SPL I T= , F6 . 3 , 3X , 1 1 HEFF I 
1C I ENCY = » F6. 3, 3X» 9HRE ACT! GN = » F6. 3, / ) 

59 FORMAT (1H0.15H FIRST STAGE N=. 12. 3X, HHWORK SPL I T=, F6 . 3 , 3X , 1 1HEFF 
1ICIENCY=,F6.3,3X,9HREACT 10N=,F6.3,/) 

60 FORMAT (1H , 30HINLET DENSITY AT HUB , ME AN , T I P= , 3F 10 . 3 ) 

61 FORMAT (1H ,29HEXIT DENSITY AT HUB , ME AN , T I P= , 3F 10 . 3 ) 

62 FORMAT ( 1H , 37HEX I T STATIC PRESSURE AT HUB , MEAN, T l P= , 3F 10 . 3 ) 

63 FORMAT (1H , 38HINLET STATIC PRESSURE AT HUB , ME AN , T IP= , 3F 10. 3 ) 

64 FORMAT ( 1H , 41HINLET STATIC TEMPERATURE AT HUB .ME AN , T I P* , 3F 1C . 3 ) 

65 -FORMAT (1H ,40HEXIT STATIC TEMPERATURE AT HUB , ME AN , T I P= , 3F 10 . 3 » 

66 FORMAT (1H , 46HI NLET TOTAL RELATIVE PRESSURE AT HUB , MEAN » T I P = , 3F' 10 

1.3) 

67 FORMAT (1H , 45HEX IT TOTAL RELATIVE PRESSURE AT HUB , ME AN, T I P= , 3F l 0 . 
13) 

68 FORMAT (1H .49HINLET TOT \L RELATIVE TEMPERATURE AT HUB , ME AN, T I P= , 3 
1F10.3) 
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69 FORMAT < 1H , 48HEX I T TOTAL RELATIVE TEMPERATURE AT HUB, ME AN, TIP** 3F 
110.3) 

70 FORMAT (1H .24HINLET TOTAL TEMPERATURE* , F12 .3 ) 

71 FORMAT (1H » 21HINLET TOTAL PRESSURE* , F 12 .3.1 

72 FORMAT (1H , 23HEX I T TOTAL TEMPERATURE*, F 12. 3 > 

73 FORMAT (1H , 20HEX I T TOTAL PRESSURE* . F 12. 3 ) 

74 FORMAT ( 1H0, 3X, 4HRA0E ,6X ,3H3 A£ , 7X , 3HVAE, 7X , 5HVEVCR , 5X » 5HUE VCR , 5X , 5 
1HWEWCR,5X,4HRADI ,6X, 3HBA I, 7X , 3HVA I , 7X ,5HVI VCR, 5 X, 5 HU I VCR ,5X,5HWIWC 
2R,2X,3HVXI,3X,3HVXE) 

75 FORMAT (1H , F8 . 3, 2F10 . 3 , 3F 10 . 4, 3F 10 . 3, 3F10 . 4.2F6 . 3 ) 

76 FORMAT (1H , 28X, 3F 10 . 4 , 3CX , 3F 10.4 1 

77 FORMAT (1H0.62H VUE ITERATION. NO SOLUTION. CHANGE RHRT OR EXTEND 
1 VUE LIMITS ) 

78 FORMAT 1 1 HO , 47HE X IT TIP DIAMETER MUST BE CHANGED. NO SOLUTION.) 

79 FORMAT 1 1H0 » 26H AREA RAT IO GREATER THAN 1) 

END 


$IBFTC ROO 

SUBROUTINE ROOT ( WH2 , KH2 1, RHRT l » RHRT2, WC , RHRT , TT, P T,GAM, CP , WF ,R, VU 
1,V,RT,RM,RH,RAPS»CI»VX) 

DIMENSION VU ( 20 ) , V(20), U C 20 ) 

G*32.17 

HA=RHRT2-RHRT1 

HAN=HA 

RALT *RHRT 1 

1 AB=WC-WH2 
C*WH21-WH2 

IF ( ABS (HAN) .LT.0.000C01 ) HAN*0.0 
H* AB/C*HAN 
RHR T=RALT + H 

A=3.1416*4.0*RM**2*( 1 .O-RHRT ) / ( 1 .O+RHRT ) 

RT=SQRT ( A/ ( 3 . 1416* ( 1 .C-RHRT**2 ) ) ) 

rh*rt*rhrt 

VUI l)*CI/RT 
VU ( 2 ) *C I /RM 
VU ( 3 ) =C I / RH 
GAMMA*GAM 

CALL DENR ( VX,G, GAMMA , CP ,R,TT,Pt»VUl2) , NF , A, RHOM ) 

CALL VXR ( VU, V, VX, WF,TT,PT,G, GAMMA, CP, R,RT,RM,RH> 

U ( 3 ) =RH*RAPS 

WHI2*( VUI3J-UI3) ) **2+VX**2 

IF ( ABS ( ( WH2 1-WH2) *. 00 1 ) .LE.C.09) GO TO 2 

IF ( ABS { WH I 2-WC ) • LE. ABS ( WH21-WH2 ) * .001 ) RETURN 

GO TO 3 

2 IF (ABS (WHI2-WC) .LE.0.09) RETURN 

3 IF ( AB * ( WH 1 2- WC ) , L T . 0 . 0 ) GO TO 4 
HAN=H 

WH2 1=WH I 2 
GO TO 1 

4 HAN=HAN-H 
WH2=WH I 2 
RAL T=RHRT 
GO TO 1 
END 
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$ I BFTC RHRT 


SUBROUTINE RHRT2 ( RHRT 1 1 , GAMMA , CP , G, R, T TI , PT I , WF , RM I , C I , * ) 
COMMON VXIL 
DIMENSION VUII3) 

GO TO 2 

1 RHRTll=RHRTI 1-.02 

IF (RHRTI 1.EQ.0. I VX I L=VXI L+ . 01 
IF ( RHRTI l.EQ.O. I RETURN 1 

2 A I = 3. 1 4 1 6*4. 0*RM 1**2* ( 1 . 0-RHRT 1 1 } / ( 1 . 0+RHRTI I ) 

RTI -SORT! At / C 3. 1416* C l .O-RHRT 1 1**2 ) ) ) 

RH I =RT I *RHR III 
VUIUMCI/RT1 
VUI (2) =CI/RMI 
VUI t 3) -CI/RHI 

CALL DENR I VX I ,G ,GAMM A »C P, R , TT I , PT l , VU I (2) ,WF, Al.RHOM) 

IF (VXI .GT. VXIL* SORT (2.0*GAMMA*R*G*TTI/(GAMMA+1.0) ) ) GO TO 1 

RETURN 

END 


$IBFTC DEN 

SUBROUTINE DENR < VX, G , GAM, CP , R , TT, PT , VUM , WF » A , RHOM ) 
VX=0.7*S0RT(2.0*GAM*R*G*TT/ (GAM+l.O) ) 

VX1=VX 
C J=778. 

1 VM=VUM**2+VX**2 

T= T T-VM / (2.0*G*CJ*CP ) 

IF IT.LT.O.u) RETURN 
P=PT*(T/TT)**(GAM/(GAM-1.0) ) 

RH0M=P/(R*T> 

VX=WF/(RHOM*A) 

IF (ABS(VX/VXl-l.O).LT. 0.001 ) RETURN 

VX1=VX 

GO TO 1 

END 


$ IBFTC VX 

SUBROUTINE VXR { VU , V , VX, wF , TT , PT , G , GAM, C P , R , RT E , RME , RHE ) 

DIMENSION VU ( 20 ) , V ( 20 ) * T(20>, P(20), RH0(20) 

VX1=VX 

CJ=778. 

1 CO 2 1=1,3 

V< I ) = VX**2 +VU ( I > **2 
T ( I )=TT-V ( I ) /(2.0*G*CJ*CP) 

P< I )=PT*( TII )/TT )**( GAM/ (GAM- 1.0) ) 

2 RHO! I ) =P ( 1 ) / ( R*T ( I ) ) 

F AB= ( RHO( 3 )-RH0 ( 2 ) )/ (RHE-RME) 

FBC=(RH0(2)-RH0(1) ) / ( RME-RTE ) 

FABC=(FAB-FBC)/(RHE-RTE) 

HTINT=RH0(3)*(RTE**2-RHE**2) /2.0+FAB*( ( R TE*» 3-RHE**3 ) / 3. 0- l RHE*RTE 



1«*2-RHE**3)/2.C)+FABC*( ( RT E*« A-RHE**A ) /A. - ( RME+RHE ) * ( RTE**3-RHE* *3 
2)/3.0+RME*RHE*(RTE**2-RHE**2 )/2.0) 

VX=WF/ (2.0*3.1A16*HTINT) 

IF {ABS(VX/VX1-1. 01. LT. 0.0005) RETURN 

VX1=VX 

GO 10 I 

END 


$IBFTC MAXI 

SUBROUTINE MAXIM ( WA , WAM AX, WTFL » I NO , TOLER ) 

DIMENSION S PEED ( 3 ) » WE IGHT ( 3 ) 

GO TO ( I, 3,5, 1A, 23), IND 

1 JUMP=1 

IF (WAMAX.LT. WA) GO TO 2 

SPEED! I ) = WA 

WE IGHT ( 1 ) =WTFL 

WA=WAMAX 

I ND = 2 

RETURN 

2 SPEED( 3 ) =WA 
WEIGHT! 3 ) = W IFL 
WA=wAMAX 

IND = 5 
RETURN 

3 SPEED! 3 )=WA 

WE I GHT ( 3 ) =WTFL 

A WA=( SPEED! I) +SPEED(3 ) )/2. 

IF (ISPEED(3)~SPEED( I)). LT .TOLER) GO TO 2A 

IND = 3 

RETURN 

5 SPEED! 2 )=W4 
WEIGHT! 2)=WTFL 

IF (WTFL. LF.WEIGHT(l) . OR. WTFL. LE. WEIGHT ( 3 ) ) GO TO 12 
I ND=A 

6 IF (WEI GHT ! 3 ) -WE IGHT ( 1 ) ) 8,9,7 

7 WAMSPEED ( 1 ) +SPEEDI2 ) ) /2.0 
RETURN 

8 WA=(SPEED!3)+SPEED(2) )/2.0 
RETURN 

9 GO TO (10, ID, JUMP 

10 JUMP=2 
GO TO 7 

11 JUMP=1 
GO 10 8 

12 IF (WEIGHT(3) .GT. WEIGHT! 1) ) GO TO 13 
WEI GHT ( 3 ) =W TFL 

SPEED! 3 ) =W A 
GO TO A 

13 WEIGHT ( 1 )®WTFL 
SPEED! 1 )=WA 

GO TO A 

1A IF ( (SPEED(3)-SPEED!1)).LT. TOLER) GO TO 2A 
IF (WTFL-WEIGHT (2) ) 18,21,15 
15 IF ( WA-SPEFD ( 2 ) ) 17,16,16 
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16 SPEED! 1 ) = SP&ED<2) 

SPEED! 2 ) = WA 
WEIGHT (1)= WEIGHT (2) 

WEIGHT I 2)*WTFL 

GO TO 6 

17 SPEED! 3 ) = SPEED(2 ) 

SPEECH 2 )=WA 
WEIGHT(3)=WEIGHT(2) 

WEI GHT { 2 ) =WTFL 

GO 10 6 

IB IF ( WA-SPEED (2 ) ) 20,19,19 

19 WE I GHT ( 3 ) =WTFL 
SPEED! 3 ) =W A 

GO TO 6 

20 WE I GHT ( l ) =WTFL 
SPEED! 1 !=WA 

GO TO 6 

21 IF ! WA . GT . SPEED ! 2 ) ) GC TO 22 
SPEEDI3)=SP£ED(2) 

WEI GHT (3)=WEIGHT(2) 

GO TO 5 

22 SPEED! 1 ) = SPE£D ( 2 ) 

WEIGHT l 1) =WE IGHT (2) 

GO TO 5 

23 SPEtDm=WA 
WEIGHT 1 1 ) = WFFL 
GO TO 4 

24 I ND=6 
RETURN 
END 


SIBFTC SUR 

SUBROUTINE SURC 

COMMON /SV/ XV A I { 5 ) , X B A I (5 ) , XBAE ! 5) , XVAE I 5 ) , SAVE ,BHS ( 5 ) , BHR ( 5 ) ,RAE 
1.(20) , TSAI 20) .WF.RAT! 3 ),RAH! 3) ,RHOI ( 10) ,RHOE( 10) , VEL ! 10 ) , V I LI 10 ) , PP 
2!10),TT(10) , C0NST2, CONST 1. V I LC 1(10), GAMMA, R, DTK i,DTR2, FAR, WETO.VXl 
32,VXE2,VXI1,VXE1 

COMMON /BVE/ BV!10,12),AB(2),LB(2),X!2) 

COMMON /CF/ BCS! 10 ) , BCR ( 10 ) , TBS AS! 10),TBSAR(10), XMRC 
DIMENSION PS I SI 15), DA SI (5 ) , PSIS2I5), DAS2(5), PSIR1 ( 5) * DAR 1(5), 
1 PSIR2I5), DAR2I5), D ALP AR ( 5 ) , DALPAS(5» 

DIMENSION B X ( 10,2) 

DIMENSION H2M(1C) , PSIM(IO), ALPHA(IO), RMCLCX(IO), B1SI5), B2SI5) 
1, B 1 R 1 5 ) » B 2 R ( 5 ) » ALPHZS (5 ) , ALPHZR(5), PS I MS ( 5 ) , PS I MR 1 5 ) * ALPHAS 
2(5), AL PHAR ( 5 ) , XCLCXS15), XCLCXR ( 5 ) » ARS(5), ARR ( 5 ) , CXS(5), CXR( 
35), 8S AS ( 5 ) , BSAR ( 5 ) , PS(5), PR!5), BPRS(5), BPRRL5) 

DIMENSION A ( 55) , B( 89 ) 

EQUIVALENCE ( BX , BV ( 1 , 1 1 ) ) 

DATA A/ 30. ,60., 90. ,120., 150., 160. , 10. , 14. , 20 . , 30. , 45. , 60 . , 70. ,79. , 

175.. 69. .60. .48. 5. 41. .37. 5. 73. 5. 69. 5. 64. .55. .43. .30. .22. .70. 5. 64. .5 

26. 5. 46. 5. 34.. 22. .14. .7?., 60., 48. ,35. ,20. ,6. 5,— 1., 82. ,52. ,35. 5, 20., 

33. , -11. ,-19. ,90., 48. ,29.5, 14., -3., -17. 5,-25. 5/ 

DATA B/0. ,20, ,40. ,60. ,80., 100. , 120. , 140. ,-20. ,-10. ,0.0, 10. ,20. ,30. 

1.40.. 50. » 60. , 1.06, 1.015, 1.0, 1.02, 1 . 075, 1 . 1 55, l .305, 1 . 57, 1. 985 , 1 . 07 
2,1.03, 1.02,1.0 3, 1.075, 1. 15 5, 1.305, 1.57,2.00,1.11, 1.07,1.06, 1.0 7, 1. 
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311, L. 13,1. M 5, 1.60, 2. 1,1.16,1.10,1.09,1.105,1.15,1.225,1.38,1.68,2 
6.16,1.18, 1.16, 1. 13,1. 16, 1.20,1.295, 1.67, 1. 76,2.18,1.26,1. 19,1.18,1 
5.20,1.26, 1.37, 1. 56 , l. . 80, 2. 22, 1. 38. 1..3C, 1.2 7, 1.29, 1.36, 1.68, 1.66,1. 
693,2.30, 1.6*. , 1.62, 1.63, 1.62,1.50, 1.61, 1.80,2.08,2.66/ 

RE in I NO 7 
1 TOT 1=5 5 
1 = 6 
J= 7 

WRITE! 7) JACK') ,K=l, I) 

I J = D J 
IP=l+l 

WRITE! 7) ( A { K ) ,K=IP,I J ) 

KK= I+J+ 1 

WRI TEl 7 M A t K ) » K=KK , I T0T1 ) 

IT0T2=89 
1 = 8 
J = 9 

WRI r E C 7 ) (B(K) , K= 1 , 1 ) 

IJ=I+J 

IP=I+1 

WRITE17) ( BIK) ,K=IP, I J ) 

KK= I +J+ 1 

WRITFI7) C B C K ) , K = KK , ITCT2 ) 

REWIND 7 
WRITE 16,19) 

VA I L=XVAI l 1 ) 

E A I L = XBA I ID 
BAEI.=XB AE I 1 ) 

VAfc I =XVAE l 2 ) 

BAEI=XBAE (?) 

BA I I=XBAI 12 ) 

VA I I=XVAI (2) 

C COMPUTE BLADE ENTERING AND LEAVING ANGLES 

B1 S ( 1 ) =90 . 

82S 1 1 ) =90 . -V A I I 
B1R ( 1 ) =90 . +BA I I 
B2RI 1 ) =90 . +B Afc I 
B1S(2)=90.-VAEI 
B2SI2) =90. -VAIL 
B1R(2)=90.+BAIL 
B2R(2)=90.+BAEL 
C 8 1 = B 1 $ 1 1 1-02 5 1 1 ) 

DB2=81R( 1 ) -82R I 1 ) 

DB3-B1S 1 2 ) -8 2S I 2 ) 

086=B1RI2)-B2R(2) 

C READ IN DATA FOR FIGURE 2 

1=6 
J = 7 

LB( 1 )=I 
L B ( 2 ) = J 

READ (7) ( B X ( K * 1) » K * 1 » I ) 

READ (7) (BXCK.2) ,K=1,J) 

READ (7) ( ( BV ( K ,L ) , L= 1 , J ) , K= 1 , I ) 

ALPHZSI 1)=BVTL(B1SII ) ,B2S( 1) ) 

ALPHZRI 1)=BVTL (B1R(1 ) ,B2R( 1)1 
ALPHZSI 2) =BV TUB IS (2) ,B2SI 2) ) 

ALPHZR(2)=BVTL(B1R(2) ,B2i<(2) ) 

C READ IN DATA FPR FIGURE 3 

NPTSB=2 
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82M l 1) =38 . 6 
PS I M I 1 ) =1 . 1 
B2M 12) = 10 . 

PSIM(2)=.26 
1 CONTINUE 

C RE AO IN DATA FOR FIGURE 5 

1 = 8 
J=7 

LB ( l ) = I 
L B ( 2 1= J 

READ 17) ( BX ( K» II * K=1 , 1 ) 

READ (7) (BX(K,2),K=1,J) 

READ (7) ((BV(K,L),L=l,J),K=l,n 

CALL LINT ( B2 S ( 1 ) » NP T SB , B2M , PS I M » PS IMS ID) 

PS I MR! 1) = I. 

CALL LINT (B2SI2) , NPTS B , B2M , PS I M , PS I MS ( 2 ) ) 
PS I 2 > = 1 . 

C READ IN DATA FOR FIGURE 4A 

PSIS1I 1 )=.4C 
PSISlt 21*1. 

DAS 1(1) =-33. 

DASI(2)=10. 

PSIS2! 1 )=.40 
PS I S2I 2 ) *1 . 

DAS 2 ! 1 I =-18. 5 
DAS2(2)=-2.5 

C READ IN DATA FOR FIGURE 48 

PSIRII l ) = . 6 0 
PS IR1I 2 ) =L . 5 
DARI 1 1 1 =—14. 5 
DARK2) -10.5 
P S I R 2 ( 1 ) = • 7 0 
PSIR2I 2 > = 1.1 
DAR 2 1 1 1 =— 9. 

0AR2 (2 ) =-6. 

CALL LINT t PS IMS ( 1 ) » 2»PS IS 1 • DAS1 »DALPAS I X) ) 
CALL LINT ( PS IMS ( 2 1 * 2 « PS I S2 » DAS2* DALPAS ( 2 ) ) 
CALL LINT (PSI MR { 1 ) » 2 »PS IR1 » DARI, DAL PARI 1 ) ) 
CALL LINT ( PS I MR ( 2 ) » 2 »PS IR2 » DAR2, DALP AR ( 2 1 ) 
DO 6 1=1,2 

6 ALPHAS! I )=ALPHZSU) + DALPAS(I) 

DO 7 1=1,2 

7 ALPHAR! I)=ALPHZR( D + DALPARd) 

XCLCXS! 1)*BVTL(DB1, ALPHAS! 1) ) 

XCLCXR ! 11=BVTL(DB2. ALPHAR! I) ) 

XCLCXS! 2) =BVTL(DB3, ALPHAS! 2) ) 

XCLCXR ( 2 )=BVTL (DBA, ALPHAR! 2) ) 

XCLCXS! l)=1.056*XCLCXSIl) 

ARS ( 1 ) =2. 

ARR! 11=3. 

A R S ( 2 1= 3 • 

ARR ! 2 ) =5 . 

FAC=1. 

8 CXS ! 1 ) =BHS ( 1 ) / ( ARS (1) »FAC 1 
IF (CXSll).LT.l.) FAC=FAC-.l 
IF (CXS(l).LT.l. ) GO TO 8 
CXR ! 1 ) =BHR ( 1 ) / (ARR ( 1 ) *FAC) 

IF (CXR(ll.LT.l. ) FAC = FAC-.l 
IF (CXR! 1 ) .LT. 1. ) GO TO t 
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CXS(2)=BHS(2 I / ( ARS (2 ) *FAC) 

IF (CXS (2) .LT. 1. ) FAC=FAC-.l 
IF (CXS ( 2 ) . LT. 1. ) GO TO 8 
CXR ( 2) - BHR ( 2 ) / ( ARR ( 2 ) *FAC) 

IF (CXR ( 2 ) .LT. 1 . ) FAC=FAC- . I 
IF ( CXR ( 2 ) . L T . 1 . ) GO TO 8 
DO 9 1=1,2 

9 BCS( I )-XCLCXS( I)*CXS( I ) 

DO 10 1=1,2 

10 BCR (I ) =XCLCXR ( I ) *CXR { I ) 

DO 11 1=1,2 

11 BSAS ( 1 ) -BHS( I ) *XCLCXS ( I ) *CXS (I ) *2.030 
DO 12 1=1.2 

12 BSARII ) =BHR ( I ) *XCLCXR 1 I ) *CXR II) *2. 030 
CO 13 1=1,2 

PS( I> = PS!MS( I)*CXS(I l/(2.*SIN(.0175*82S( I ) ) *SIN( .0 175*B2S ( I ) > * ( COT 
1 AN( ,0175*B2S( I ) ) -COT AN ( . Cl 75*B1S ( 1)1)) 

13 CONTINUE 

DO 14 1=1,2 

PR (I ) -PS I MR ( I )*CXR(I ) / 1 2 .* S I N ( . 01 75*B2R ( I ) )*SIN( ,0175*B2R( I) )*(COT 
1 AN ( .0175*B2R( I ) ) -COT AN ( . 01 75 *B 1R ( I ) ) ) ) 

14 CONTINUE 

XMRC=2. *3. 14159*RAE( 2) 

DO 15 1=1,2 

15 BPRS ( I )= XMRC/PS ( I ) 

DO 16 1=1,2 

16 BPRRII )=XMRC/PR{ I ) 

DO 17 1=1,2 

17 TBS AS ( I )-BPRS ( I ) *BSAS ( I ) 

DO 18 1=1,2 

18 TBSARt I )=BPRR( I )*BSAR( I) 

WRITE (6,20) BIS ( 1 ) , B1R ( 1 ) * B IS ( 2 ) »B1R ( 2 ) 

WRITE (6,21) B2S(1),B2R(1),82S(2),B2R(2) 

WRITE (6,22) ALPHAS ( 1 ) , ALPHAR ( 1 ) » ALPHAS ( 2 ) , ALPHAR { 2 ) 

WRITE (6,23) XCLCXSI 1 ) , XCLCXR ( 1 ) , XCLCXS ( 2 ) , XCLCXR ( 2 ) 

WRITE (6,24) CXS (1)»CXR( 1 ) , CXS ( 2 ) , CXR ( 2 ) 

WRITE (6,25) BSAS(l) , 8SAR( 1 ) , BSAS (2) ,BSAR(2) 

WRITE (6,26) PS ( 1 ) ,PR ( 1 ) ,PS(2).PR(2) 

WRITE (6,27) XMRC 

WRITE (6,23) BPRS ( 1 ) , BPRR( 1 ),» BPRS ( 2 ) , BPRR ( 2 ) 

WRITE (6,29) TBSAS ( 1 ) ,TB SAR ( 1 ) , TBS AS ( 2) » TBSARt 2 ) 

CALL COOLA 
RETURN 
C 
C 

19 FORMAT (1H1.54HINITI AL SURFACE AREA CALCULATIONS ***************** 
1***,//) 

20 FORMAT (1H » 45HBL ADE ENTERING ANGLE, BETA l.FOR ROWS 1 ,2 , 3 , 4= , 4E 1 2 . 
13) 

21 FORMAT (1H , 44HBL ACE LEAVING ANGLE, BETA 2,F0R ROWS 1 , 2 , 3 , 4= , 4F 12 . 3 
1) 

22 FORMAT (1H , 40H6L ADE CAMBER ANGLE , PH I , FOR ROWS 1 , 2 , 3 , 4=, 4F 12. 3 ) 

23 FORMAT (1H , 69HBL ADE CAMBER LENGTH TO CHORD LENGTH RATIO FOR ROWS 

11,2,3, 4=»4F12. 3) 

24 FORMAT (1H , 30HAX I AL CHORDS FOR ROWS 1 , 2 , 3 ,4=, 4F12 . 3 ) 

25 FORMAT (1H , 36HSINGLE BLADE AREAS FOR ROWS 1 ,2,3,4=,4F12.3) 

26 FORMAT (1H » 23HP ITCH FOR ROWS 1 , 2 , 3, 4=, 4F 12 .3 ) 

27 FORMAT (1H , 26HMEAN RADIAL C IRCUMFfcRENCE= , F 12.3) 

28 FORMAT (1H , 33HNUMBE R OF BLADES IN ROWS 1, 2 , 3, 4= , 4F 12 . 0 ) 
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ui o o ^ uj rv) i— no non 


29 FORMAT ( IH » 42HT0TAL BLADE SURFACE AREAS IN ROWS 1 , 2 , 3, 4= , 4F 1 2 . 3 ) 

30 FORMAT 1215) 

31 FORMAT I8F10.1J 

32 FORMAT (15) 

33 FORMAT (2F12.4) 

END 


$ I B F T C BVD 

FUNCTION BVTL (V1,V2) 

COMMON /BVE/ BVI 10,12) ,AB( 2 ) * LB( 2 ) *X ( 2 ) 

DIMENSION BX (10,2) 

EQUIVALENCE ( BX , BV ( 1 , 1 1 ) ) 

BIVARIATE INTERPOLATION FOUR POINT FORMULA 

X ( 1 >=V1 
X l 2 ) =V2 
DO 4 1=1,2 

IS ARG BELOW RANGE OF TABLE 
IF (XH ).LT.BX(1,I)) GO TO 2 
NO 

N=L 8 ( I ) 

DO 1 J=2,N 

IF (X(I).LE.BX(J,D) GO TO 3 
WR I TE (6,5) X ( 1 ) , X ( 2 ) 

CALL EXIT 
M= J-l 

IF (I.EQ.l) K=M 

AB ( I ) = ( X ( I ) — B X ( M t I )) /(BX (J, I >-BX(M, I ) ) 

I =K + 1 

BVTL=BV(K,M)+(BV(I,M)-BV(K,M) ) * AB ( 1 ) + ( BV ( K , J ) -BV ( K, M ) ) * AB { 2 ) + ( BV ( K 
i,M)+BV( I, JI-BVl I,M)-BV(K,J) ) *AB( 1 )»AB(2) 

RETURN 


FORMAT ( 1 H 1 , 27 HB I V AR I A TE INTERP ERROR Vl= ,610.2, 5H V2= ,E10.2) 
END 


$ IBFTC OlfNT 

SUBROUTINE LINT ( XX, NT AB , XT AB , YT AB , Y Y ) 
DIMENSION XTAB(l), YTAB(I) 

DO 2 N = 2» NTA6 

IF ( XX-XTAB ( N ) ) 1,1,2 

1 N2 = N 
N1=N-1 
GO TO 3 

2 CONTINUE 
N2=NTAB 
N1=N2-1 

3 IF (YTAB(N1)~YTABIN2) ) 5,4,5 
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4 V Y= YT&B (Ml) 

RETURN 

5 IF (XTAB( N1)-XTAB(N2 ) 1 7,6,7 

6 YY=0.5» (YTA6(N1)+YTAB(N2 J) 

RETURN 

7 B=(YTAB(N2)-YTA8(NL) ) / ( XTAB ( N2 ) -XTAB ( NI ) I 
A=YTAB(N1 )-B*xtAB(Nl ) 

YY= A+B*XX 
RETURN 

END 


tIBFTC COOL DEBUG 


1 


SUBROUTINE COOLA 
REAL K 

CORDON /X LOUP/ RY,RYB,RWC,RDHT,NWC, ITT, TESTY 

COMMON /CFCT/ TCI , CORRFW.CPRRF, YW,TPPMTP( 100) , UHC , WPS TAT < 5 > ,HTtMP( 
15) 

COMMON /SV/ XVA I (5), XBAI 1 5 ) , XB A£ ( 5 » , XV AE ( 5 ) , S A VE , BH S ( 5 ) , BHR ( 5 > , R AE 
1(20) , TSAI 20) ,WF,RAT( 3) ,RAH(3 ) ,RHGI ( 10) ,RHOE( 10) ,VELI 10 ) , V I L ( 10 ) , PP 
2(10) ,TT(10) , CUNST2, CONST 1, VILCI 1 10 ) , GAMMA, R, DTR1 ,DTR2 , F AR , WE TO , VX I 
32*VXE2,VXI1,VXE1 

COMMON /CF/ BCSl 10), BCR( 10) , TBSASI 10) , TBSAR( 10) ,XMRC 
DIMENSION RC ( 20 ) » BS A ( 20 ) , TG(20), R I ( 2 C ) , RE(20), RVA(20), X(20), 
1 V I SC ( 20 ) , RE A8 ( 20 ) , CONDI 20), HG(20), U(20), TMS(20), TMR(20), TM 
2A ( 2 0,20 ) * OR ( 2'J * 20 ) » T.M.w 1CM ( 10,10) , WCWE(10,10) , TPN(IO), TWF(IO), 
3 TWrO(lO), BSAO(IO), RVO(IO), WL(IO), VCR(IO), VOVCR(IO), R.ORP(IO) 
4, RV(10), VX I ( 4 ) , VXE ( 4 ) , TPRI20), PR(20), PRN(20), TCON(20), CON( 
520), TMU120), V I S ( 20 ) , T DARK 1C ) , EF(lu), CXF(lO), CPG(IO), RATCP( 
610), BRAKT(IO), WCWEC(IO), TGARTAC1C) , PRW(IO), CONDG(IQ), PRG(IC) 
7, V I SC S ( 9 ) , TPPN ( 50 ) , TGAMAX(IO) 

DIMENSION BH ( 5 ) » WBH(5), WCW(5), wRHO ( 5 ) , WVEX05), WCWT<5), WDA(5) 
1, WANNA { 5 ) , WBSA ( 5 ) 

DIMENSION WCR(IO) 

DIMENSION Tfc STC ( 50 ) 

DIMENSION TTWF(IO) 

IF( ITT.GT.DGO TO 5001 

READ(5,43)FARC,PF,DTY 

PTD=0. 

READ (5,42) NP TO 
NPTPR=NPTD 


NPTMU=NPTD 
NPTCON = NP TO 


CO 1 1 = 1 , N P T D 

READ (5,32) TPR ( I ) ,PRN ( I ) , V I S ( I ) »CON ( [ ) 
TMU ( I) =TPR ( I ) 

TCOMI I = T P R { I ) 

CONTINUE 
NTS = 1 

READ (5,43) (TMS( I), 1=1, NTS) 

NTR = 1 


READ (5,43) 
READ (5,43) 
READ (5,43) 
READ (5,52) 


(THRU), 1 = 1, MR) 

THERME 

CPC 

TwALLA 
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O CD 


5001 

2 


3 

4 

5 

6 

7 

9 

10 


TP2=TC I 
NSTAGE=2 
WRI TE( 6,41 > 

WRITE (6,33) 

WRITE (6,34) PTD»FARC»PF,DTY 
WRITE (6,35) NPTD 
00 2 1=1, NPTD 

WRITE (6,36) TPR( I ), PRN( !) ,VIS( I ) ,CGMC I) 

WRITE (6,37) NTS, TMS ( 1 ) * NTR » TMR ( l ) 

WRITE (6,38) THERMEt CPC, TWALLA 
NR0WS=2*N$TAGE 
NR0WSP=NR0WS+1 
G = 3 2 . 2 

BC ( 1 ) =RCS ( l ) 

BC ( 2 ) = BCR ( l) 

BC ( 3 ) = B C S ( 2 ) 

BC ( 4 ) = BCR ( 2 ) 

6 S A ( 1 ) = T B S AS ( 1 ) 

8SA(2)=TBSAR( 1) 

BS A ( 3 ) = TB S A S ( 2 ) 

BSA(4)=TBSAR(2) 

WRITE (6,47) 3 S A ( 1 ) , B S A ( 2 ) ,BSA( 3) , B S A ( 4 ) 

WBSA(l)=BSA( 1) 

W8SM2) =3SA(2) 

W8S A ( 3 ) =8SA( 3 ) 

WBSA(4)=BSA(4> 

R I ( 1 )=WF*144.*C0RRFW/ (3. 14 159* ( RAT ( i ) **2 . -R AH ( i ) ** 2 . ) ) 
AX= 3 • 14 159* ( RAT ( 1 ) ** 2 . -R AH ( 1 ) **2 • ) 

CO 3 1=2, NRUWS 

R I ( I ) = R HO I ( I ) * V I L ( I) *CORRF 

CONTINUE 

DO 4 1=1, NROWS 

RE ( I )=RHOE( I ) *VEL ( I) *CORRF 

CONTINUE 

TPN ( l ) = TGA ( 1 ) 

TWF ( 1) =WF*CORRFW 
CO 5 1=1, NROWS 
WC( I 1 = 3.0 
DO 6 1=2, NROWSP 
TWF 1 1) =TWF ( l-l ) 

TPN ( I ) = TT ( I ) 

WE=WETJ*COR«FW 
WRITE (6,44) TCI 
WRITE (6,39) THERME 

WRITE (6,49) HTEMPU I ,HTEMP(2) ,HTEHP(3) .HTEMP14) 

WRITE (6,46) TGA ( 1 ), TGA( 2), TGA (3), TGA (4) 

CONTINUE 

TGAMAXU ) =TGA( 1 ) +PF* ( TGA ( 1 ) - ( TP2+DT Y ) ) 

TGA?' 1 AX ( 3 ) = TO A ( 3)+(TGAMAX(l ) - TGA ( 1 ) )/2. 

CO 3 1=1, NROWS 

RVA ( I)=(RI ( I ) +RE ( I ) ) / 2. 

CONTINUE 

STORE VALUES FROM PREVIOUS ITERATION 
DO 9 1=1, NROWSP 
TWFC ( I ) = T W F ( I ) 

CONTINUE 

CO 10 1=1, NROWS 

BSAO( I ) =B S A { I ) 

RVO ( I ) — ( R I (I ) +RE ( I 0/2. 
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AXO=AX 

DO 11 1=1 , NROWS 
XIII =BC C I ) / 1 2 . 

11 CONTINUE 

DO 12 1=1 » NROWS 

CALL LINT (TGA(I) , NPTMU , TMU , V I S , V I SC ( 1)1 
REA8II ) = ( RVA ( I)*X( I) /VISCI !>)**. 8 

12 CONTINUE 

00 13 1=1 tNROWS 

CALL LINT ( TGA ( I ), NPTCON ,TCON, CON, CONDI I I ) 

CALL LINT ( TGA( I) , NPT PR ,TPR, PRN.PR (I) ) 

C COMMUTE BLADE OUTSIDE HEAT TRANSFER COEFFICIENT 

HGl I )=• 037* I PR 1 1) **• 33 ) *COND ( I ) *REA8 ( I ) /XI I ) 

13 CONTINUE 

DO 14 I = 1 » NROWS 

U( I I =HS ( I I*BSA( U /144. 

14 CONTINUE 

CO 15 J = 1 » NTS 
TMAU, J)=TMS( JJ-PTD 
TMA ( 3, J ) =TMS I J ) -PTD 

15 CONTINUE 

DO 16 J = 1 , N I R 
TMA(2»J)=TMR(J)-PTD 
TMAI4, J)=TMR{ JJ-PTD 

16 CONTINUE 
NTM=NTS 

DO 17 1=1, NROWS 
DO 17 J = 1 » N T M 

OR ( I » J ) =U I I )*(TGA( I ) -TMA (I , J ) ) 

17 CONTINUE 

DO 18 1=1, NROWS 

DO 18 L=1 , NTM 

TMMTCMI I, L)=TMAt I,L)-TCI 

WC WE (I , L I =QR ( I , L ) / I W E * C P C* THE R M E *T M MT C M ( I , L ) ) 

IF ( WC W E ( I , L ) . L T . 0 . ) WCWCl I ,L)=0. 

WC( I )=WCWE( I ,L)*WE 

18 CONTINUE 

WC( 1) = WC1 1)*(TGAMAX( 1 )-TMA( 1 , 1 ) ) / ( TGA 1 1 l-TMA ( 1 , 1 ) ) 

WC ( 3 ) = WC ( 3 ) * ( TGAMAX ( 3 )-TMA (3 , 1 )) / I TG AT 3 ) -TMA ( 3 , 1 ) ) 

WC WE ( 1 , 1 1 = WC ( 1 ) / W E 
WC W E ( 3 » 1 ) = WC ( 3 ) / W E 
WCT=0.0 

DO 19 1=1, NROWS 

19 WCT = WCT + WC(n 
Y= ( wCT/WE ) +YW 

TWF(l)=WE*l l.-Y)*( l.+FARC) 

C AD I ABAT I CALLY MIX BLADE COOLANT FLOWS WITH GAS STRFAM 

DO 20 I =2 »NROWSP 
TWF { I ) = TWF I I -1 ) +WC ( I — 1 ) 

IF U.F0.2) TPNt I » = ( ( ( TWF{ 1-1 ) *TPN ( 1-1 ) ) )*■■( WC ( 1-1 ) *TP2 ) ) /TwF ( I ) 
TPPN (2 ) =TPN ( 2 ) +C0NST2 

IF ( I.EO. 3) TPPN ( I )=< TWFII-l )*TPPN C I — 1 ) + ( WC( I- 1 ) *TP2 ) ) /TWF(I) 
TPN ( 3 ) = TP PN { 3 ) -TPPMTP ( 3 ) 

IF ( I.EO. 4) TPN( n = ( ( (TWF( l-l ) *TPN C I — 1 ) 1 )*■( WC { I — l ) *TP2 ) ) / TWF 1 1 ) 
TPPN ( 4 )=TPN { 4 ) +CONST 1 

IF ! I.EO. 5) TPPN ( I )=( TWFM-1 )*TPPN( I - 1 1 + ( WC ( I - 1 1 * T P 2 I ) / TWF ( I ) 
TPN ( 5) = TPPN ( 5 ) -TPPMTP ( 5 ) 

VCR ( I ) =SQRT { { 2 . *GAMMA/ ( GAMMA+1. ) ) *G*R*TPN ( I ) ) 

WCR ( I ) = SORT I ( 2 . *GAMMA/ (GAMMA + 1 . ) >*G*R*TPPNM ) ) 
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VOVCR ( r ) = V ILC I ( I ) / VCR ( I ) 

RORP( I 1 = 1 1.-( (GAMMA-1. > / { GAMMA* 1. ) ) *VOVCR( I )**2. )**(!./( GAMMA- 1 . ) ) 

20 CONTINUE 

VXI ( 2) =VXI 1/ VCR ( 2 ) 

VXf. (2)=VXEl/VCR(3) 

VXI (4) =VX I 2/ VCR ( 4 ) 

VXE (4)=VXE2/VCR(5) 

RI ( I ) = TWF ( 1 )*144./AX 
00 21 I =2 » NROWS 

21 RI ( I > = ( VI L( I )*PP( I )*RORP (I )/ (R*TPN<1 ) » )*CORRF 
DO 22 1=1, NROWS 

REC I )= ( VEL ( I )*PP( 1+1 ) *RORP( 1+1 )/(R*TPN( l + l ) ) )*CORRF 

22 RV( Ii=(RI ( I)+RE( I))/2. 

TGA 11) =TPN ( 1 ) 

TGA ( 2 ) = TPPN { 2 ) 

TGA ( 3 ) =TPN { 3 ) 

TG A ( 4 ) = TPPN ( 4 ) 

DO 23 1=1, NROWS 

23 BSA ( I ) =BSAO( I ) * ( TwF( I >/TWFO( I ))*{RVO( I )/RV (I)) 

AX = AXO*BSA ( 1 ) /BSAOd ) 

BH( 1)=BHS( 1 ) *BSA ( 1 1 / WBSA ( l ) 

BH ( 2 ) =BHR l 1 ) *BSA( 2I/WBSA (2 ) 

BH(3)=BHS(2)*BSA(3)/WBSA(3) 

BH ( 4 )=BHR l 2 ) *BSA (4 )/ WBSA (4 ) 

DO 24 1=1, NROWS 
C CHECK FOR SOLUTION 

24 IF (ABSIBSAU )-BSAO( I ) ).GT.. 00001) GO TO 7 
WRITE (6,48) TGA(1),TGA(2) , TGA ( 3 ) , TGA ( 4 ) 

WRITE (6,50) TPN I 2 ) « TPPN (3 ) » TPN( 4 ) , TPPN { 5 ) 

WRITE (6,51) BS A ( 1 ) , B S A ( 2) , BSA ( 3 ) , BSA ( 4 ) 

DO 25 1=1, NROWS 

DO 25 L=1 ,NTM 

WRITE (6,45) I , TMA (I , L ) , WC WE ( I » L ) 

25 CONTINUE 

9000 FORMAT ( 1H .44HTURBINE WEIGHT FLOW AT EXIT OF ROWS 1 , 2 , 3, 4= , 4F12 . 6 ) 
IF ( NWC. EQ. 1 ) GO TO 300 
C COMPUTE WALL COOLANT FLOWS 

DO 26 1=1, NROWS 
TGART( I )=TGA( I )-150. 

EF ( I )= ( TWALLA-TGART ( I ) ) / ( TC I-TGART ( I ) ) 

CXF ( I ) =1 . 2*X ( I ) 

D=XMRC/(12.*3.14159) 

CALL LINT (TGARTt I ) , NPTCON , TCON ,CON , CONDG (I) ) 

CALL LINT (TGARTl I ), NPTPR, TPR,PRN, PRG( It) 

CALL LINT (TGARTU ) , NPTMU, TMU, V I S , V I SCG ( I ) ) 

CPG( I) =PRG( I ) *CQNDG( I )/V ISCGU ) 

RATCP( I ) =CPG ( I ) /CPC 
TGARTA( I ) = ( TC I +TGART ( I ) ) /2. 

CALL LINT ( TGARTA( I ) , NPT PR , T PR , PRN, PRW ( I » ) 

BRAKT ( I )=( ( 1.9*(PRW( I ) ** .667 ) ) / EF ( I ) )-l. 

VI SGG= V I SCG ( I ) 

WCWECt I ) = ( . 329*R ATCP ( I ) * ( V I SGG**. 2 ) * ( RV ( I ) * * . 8 ) * ( .8*3. 14159*D/WE ) * 
l-( CXF { I ) ** . 8 ) ) / BRAKT( I ) 

WCWEC ( I ) =2 . *WCWEC ( I) 

IF (WCWEC II) .LT.O.O) WCWEC(I)=0.0 
WCW ( I ) =WCWEC ( I ) *WE 

26 CONTINUE 

WRITE (6,40 TWALLA 
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WCWECT = 0» 

DO 27 1=1* NROIrt S 
WCWECT=WCWECT+WCWEC( I ) 

27 CONTINUE 

DO 28 1 = 1 » NROWS 
WRITE (6,53) I,WCWEC!I) 

28 CONTINUE 
RTCI =R*TC I 

WRHOM ) =WPST ATI 11/ RT Cl 
WRHOI2) =WPSTAT(2)/RTCI 
WRHQ(3)=WPSTAT(3)/RTCI 
WRH0(4)=WPSTAT(4)/RTCI 

wvexm=vxn/ 2 , 

WVEX(2)=VXEl/2. 

WVEX(3)=VXI2/2. 

WVEX(4)=VXE2/2. 

WCWT { 1 ) =WCW ( 1 ) 

DO 29 1=2,4 

29 WCWTII )=WCWT( I-1)*WCW( 1) 

DO 30 1=1,4 

WDA(I)=!WCWT(I )/!WRHO( I>*WVEX(1) ) 1*144. 

WANNA! I )=XMRC*BH( I ) 

30 WBH (I) =BH ( 1 ) * { WANNA! I ) + WDA 1 I ) I /WANNA! I ) 

300 CONTINUE 

IF(NWC.EQ.1)WCW!1)=0.0 
IF{NWC.EQ.1)WCW!2)=0.0 
IF(NWC.EQ.l)kCW(3)=0.0 
IF ! NWC. EQ • 1 ) WCW ! 4 ) =0. 0 
TTWF(2)=TWF(2)+WCW(1 ) 

TTWF(3)=TWF(3H-kCW(l >+WCW!2) 

TTWF{4)=TWF(4)+WCW(1)+WCW(2)+WCW(3) 

TTWF (5 ) =TWF( 5 )+WCW ( 1 1 +WCk( 2 ) +WCW! 3 ) +WCW! 4 ) 

WRITE! 6,9000 ) TTWF 121, TTWF! 3 1 , TTWF I 4 » , TTWF ( 5 ) 

IF ( NWC. EQ. 1 ) WRITE ( 6, 64 ) BH! 1 ) ,BH( 2 ) » BH ( 3 ) »BH( 4) 

IF ! NWC. EQ. 2) WRITE! 6, 65) WBH! 1),WBH(2) ,WBH( 3 ) , WBH! 4) 

YB=WCT/WE 
YW=WCWECT 
IF(NWC.EQ.l) YW=0. 

Y=YW+YB 

WRITE 16,66) VXl (2)«VXE( 2) , VXI 14) , VXE14) 

IF ( NWC.EQ. 1 ) WRITE (6, 62) T PN (5 ) 

IF { NWC. EQ. 1 ) GO TO 302 

WTPF=!TPN!5)*TWF15M-TCI*WCWECT*WE)/!TWF!5)+WCWeCT*WE) 

WRITE ! 6, 63 ) WTPE 
302 CONTINUE 

WRITE ! 6 , 67 ) HG! 1 ) «HG ( 2 ) *HG ( 3) ,HG1 4) 

WRITE (6,68) VCR ( 2 ) , VCR! 3) , VCR! 4) ,VCR( 5) 

WRITE! 6,69)WCR12)fWCR(3) ,WCR!4) ,WCR(5) 

69 FORMAT ( 1H ,54HCRITICAL RELATIVE VELOCITIES AT EXITS OF ROWS 1,2,3 
1 ,4=, 4F12. 3 ) 

WRITE (6,54) 

C COMPUTE VALUES FOR RERUN IF NECESSARY 

TESTC! ITT) =ABS!Y-TESTY) 

IF(TFSTC( ITT). LE.. 00 10) WRITE (6,6000 ) 

6000 FORMAT !IH , 54HTHE LAST SET OF OUTPUT CONTAINS THE CONVERGED S0LUT1 
1CN ) 

IF (TESTC! ITT).LE..OOIO)STOP 
IF! ITT.EQ.l) GO TO 31 

IF! ABSITESTC! ITD-TESTC! ITT-1) ) .LT..0005)G0 TO 500 
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31 CONTINUE 
WON=WE*tl.+FAR*( l.-YI-Y) 

DHCB=i.+.5*l YB/( ( l.-Y)*( l. + F AR ) ) ) * ( TCI / TGA ( 1 ) ) 

DHTN=DHC/ { ( ( l.-Y)*ll.+FAR) )*OHCB> 

RYB=YB 

RY=Y 

RDHT=DHTN 

RWO=WON 

RETURN 

500 IF(Y.GT.TESTY) Y= Y-TE STC( ITTI/2, 

IF (Y.LT. TESTY) Y=Y+TESTC( ITTI/2. 

YB=Y-YW 
GO TO 31 

C 

32 FORMAT (4F12.7) 

33 FORMAT (1H , 12HINPUT DATA**) 

34 FORMAT (1H , 4HPTD= ,F8 .3, 17X, 5HFARC=,F8.3, 16X.3HPF* , F8.3, 18X,4HDTY= 

1.F8.3) 

35 FORMAT (1H ,5HNPTD=,I3) 

36 FORMAT (1H , 4HTPR= ,F8 . 3 , 17X, 4HPRN= , F8. 3 , 17X, 4HV I S= , E9 . 3, 1 7X , 4HC0N* 

1,E9.3) 

37 FORMAT ( 1H , 4HNTS= , 1 3 , 22X, 4HTMS*, F8 . 3, 17X , 4HNTR* , 1 3 , 22X, 4HTMR = , F 8 . 
13) 

38 FORMAT (1H , 7HTHERME* » F8 .3, 14X, 4HCPC*» F8 .3, 17X ,7HTWALLA= , F8. 3, /// ) 

39 FORMAT tlH ,44HASSUMED BLADE THERMAL COOLING EFFECTIVENESS=»F8.3 ) 

40 FORMAT { 1H > 31HASSUMED WALL METAL TEMPERATURE*, F12. 6) 

41 FORMAT ( 1H1 , 40HC00LANT CALCULATION ********************,//) 

42 FORMAT (13) 

43 FORMAT (6F12.6) 

44 FORMAT { 1H ,26HC00LANT INLET TEMPERATURE* , F 12 . 3 ) 

45 FORMAT (1H , 34HASSUMED METAL TEMPERATURE FOR ROW( , I2,2H)*,F12.3, 10 
1 X, 25HBL ADE COOLANT FLOW RATI0*,F12.6) 

46 FORMAT ( 1H .57HINITIAL TOTAL INLET TEMPERATURE RELATIVE TO ROWS 1, 
12»3,4*,4F12.3) 

47 FORMAT (1H ,51HINITIAL TOTAL BLADE SURFACE AREAS FOR ROWS 1,2, 3,4= 
1,4F12.3> 

48 FORMAT ( 1H , 57HREV ISED TOTAL INLET TEMPERATURE RELATIVE TO ROWS 1, 
12,3,4= ,4F12.3) 

49 FORMAT (1H ,56HINITIAL TOTAL EXIT TEMPERATURE RELATIVE TO ROWS 1,2 
1 ,3,4=, 4F12.3 ) 

50 FORMAT (1H ,56HREVISED TOTAL EXIT TEMPERATURE RELATIVE TO ROWS 1,2 
1,3,4=, 4F12. 3) 

51 FORMAT (1H ,51HREVISED TOTAL BLADE SURFACE AREAS FOR ROWS 1,2, 3, 4= 
l ,4F12. 3 ) 

52 FORMAT l E12. 3,F12.3) 

53 FORMAT C 1H , 32HWALL COOLANT FLOW RATIO FOR ROW ( , 12, 2H ) =, F 12 . 6 ) 

54 FORMAT (1H ,///, 20X, 20H* *******************,/// ) 

62 FORMAT (1H , 50HTURBI NE EXIT TOTAL TEMPERATURE WITH BLADE COOLANT*. 
1F12.3) 

63 FORMAT (1H , 59HTURBI NE EXIT TOTAL TEMPERATURE WITH BLADE AND WALL 
1C00LANT=,F12,3) 

64 FORMAT (1H , 58HAVERAGE BLADE HEIGHTS FOR ROWS 1, 2.3,4 WITH BLADE CO 
'10LANT*, 4Fl2. 3) 

65 FORMAT! 1H .67HAVERAGE BLADE HEIGHTS FOR ROWS 1,2, 3, 4 WITH BLADE AN 
ID WALL COOLANT® , 4F 12 . 3 ) 
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66 FORMAT { 1 H , 54HREVISEC V.XI AND VXE FOR FIRST AND SECOND STAGE ROTO 
1RS= , 4F12. 3 ) 

67 FORMAT <IH , 43HAVERAGE H.T. COEFFICIENTS FOR ROWS I ,2 , 3, 4* . 4E 12. 3 ) 

68 FORMAT (1H . 54HCR ITI C AL ABSOLUTE VELOCITIES AT EXITS OF ROWS 1,2,3 
1,4=,4F12.3) 

END 
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APPENDIX C 


SYMBOLS 


A 

AR 

BH 

C x 


D 

f/a 

h 

Ah 

k 

N 

P 

PF 

Pr 


P 

Qa 

Qr 

R 

RHRTI1 

RHRTI2 


r 

T 

t 

U 


area 

area ratio, A a j/ A a 2 

blade height 

blade camber length 

specific heat at constant pressure 

blade axial chord length 

diameter 

primary burner fuel-air ratio 
convective heat-transfer coefficient 
specific enthalpy change 
thermal conductivity 
turbine rotative speed 
blade pitch 

pattern factor, (T'^ max - TjJ/ff” - T^ s ) 

Prandtl number 
pressure 

available heat sink 
required heat sink 
gas constant 

FORTRAN symbol for first assumed hub-tip radius ratio 
FORTRAN symbol for second assumed hub-tip radius ratio 
turbine hub reaction 
radius 

temperature 
blade thickness 
blade speed 
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V absolute velocity 

VXEL FORTRAN symbol for limiting value on (V /V ) 

x cr 2 Iy 

VXIL FORTRAN symbol for limiting value on (V /V ) 

x cr l w 

VXEI FORTRAN symbol for . limiting value on (V /V ) 

x cr g n 

VXH FORTRAN symbol for limiting value on (V /V 0 J 

x cr 1,H 

W relative velocity 

w weight flow 

j3 angle of relative velocity from tangential 

y specific heat ratio 

77 turbine efficiency 

77 ^ blade thermal cooling effectiveness 

p density 

r stage work fraction 

(p blade camber angle 

aerodynamic (tangential) blade loading coefficient 
v constant, 77 = 3.14159 

Subscripts : 
a annulus 

BM blade metal 

b blade(s) 

c compressor 

cr critical, corresponding to conditions at Mach 1 

E engine 

g turbine main gas stream 

H value at hub radius 

i indicates substitution at HUB(H), MEAN(M), or TIP(T) 

M value at the mean radius 

max maximum value 

J-U 

n n iteration 

S surface 
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T value at tip radius 

TOT total 

t turbine 

u component in direction of blade speed vector 

w wall 

x component in turbine axial direction 

y cooling air 

Z Zweifel (from ref. 5) 

1 inlet 

2 exit 

Q first stage 

(D second stage 

I first blade row 

II second blade row 

UI third blade row 

IV fourth blade row 

Superscripts: 

average value 

” total state condition relative to rotating blading 

» total state condition relative to stationary blading 
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TABLE I. - DEFINITION OF UNITS USED IN 
COMPUTER OUTPUT 


Parameter 

Units 

Pressures 

lbf, ft 2 

Temperatures 

°R 

Velocities 

ft sec 

Densities 

lbm ft 3 

Angles 

degree 

Lengths 

inches 

Areas 

in . 2 

Convective heat- transfer coefficients 

Btu/(sec)(£t 2 )(°R) 

Turbine weight flows 

lbm/sec 
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developed turbines. Figure 7. - Velocity diagram Fortran nomenclature. 








